NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Report 32-987 

Acceleration of Liquids in 
Two-Phase Nozzles 


D. G. Elliott 

E. Weinberg 


Approved by: 


D. R. Bartz, Mcytager 

Research and Advanced Concepts Section 


JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 

July 1, 1968 


TECHN/CAI REPORT 32-987 


Copyright © 1 968 
Jet Propulsion Laboratory 
California Institute of Technology 

Prepared Under Contract No. NAS 7-100 
National Aeronautics & Space Administration 



Acknowledgment 

Important contributions were made by several people. L. G. Hays performed 
the one-component nozzle analysis. D. J. Cerini conducted the nitrogen-water 
tests on the 50-in. nozzle. H. Otte conducted the Freon-water tests. The computer 
programs were developed by C. Milton, B. Williams, A. Voskressensky, and 
P. Kaestner. 


JPL TECHNICAL REPORT 32 - 9 87 


Hi 



Contents 


I. Introduction 1 

II. Analysis 2 

A. Two-Component Flow 3 

’ 1. Assumptions 3 

2. Derivation of equations for free-stream flow 3 

3. Optimum pressure profile 10 

4. Wall shear and boundary layer 11 

| 5. Auxiliary quantities 12 

\ 6. Units 14 

! B. One-Component Flow 14 

1. Assumptions 14 

| 2. Derivation of equations for free-stream flow 14 

i 3. Optimum pressure profile 15 

4. Wall shear and boundary layer 15 

5. Auxiliary quantities 15 

6. Units 15 

C. Isentropic Flow 15 

1. Two-component flow 15 

| 2. One-component flow 16 

i 

III. Computer Programs 16 

IV. Two-Phase Nozzle Design 16 

V. Comparison with Experiment 21 

A. Nitrogen— Water Tests in the 50-in. Nozzle 21 

1. Construction and test procedure 21 

2. Pressure profile 27 

3. Flow rate 28 

4. Exit velocity 29 

B. Freon-Water Experiment 29 

1. Flowrate ..29 

- 2. Exit velocity 29 

C. Six-Inch Nozzle .32 

" 1. Flowrate 32 

2. Exit velocity 32 


IV 


JPL TECHNICAL REPORT 32-987 



Contents (contd) 


D. Steam Nozzle ^ 

E. Liquid Nitrogen Nozzle 34 

VI. Summary and Conclusions 34 

Nomenclature 35 

References 41 

Tables 

1. Conversion factors employed in computer programs 14 

2. Theoretical performance of 50-in. nozzle with water and nitrogen . . 28 

3. Measured performance of 50-in. nozzle with water and nitrogen . . 30 

4. Theoretical and measured performance of 50-in. nozzle with 

water and Freon 1301 at 150 psia 30 

Figures 

1. Two-phase nozzle flow geometry and nomenclature 2 

2. Drag coefficient of spheres and drops 9 

3. Annular nozzle geometry . 13 

4. Correspondence between nozzle length and Lagrangian 

multiplier - ^ 

5. Normalized profiles of optimum nozzles 18 

6. Modification of optimum shape to obtain experimental nozzle . . .18 

7. Theoretical pressure profile for experimental nozzle and 

tabulated values employed in real-flow analysis 19 

8. Effect of length and initial drop diameter on exit velocity, 

with and without breakup 19 

9. Variation of drop diameter with distance from inlet ..... 19 

10. Variation of drop diameter with travel time 20 

1 1 . Effect of nozzle length and flow rate on exit velocity 

with friction 20 

12. Optimum and off-optimum pressure profiles and 

corresponding exit velocities 20 

13. Variation of liquid and gas velocities with distance 21 

14. Variation of liquid and gas temperatures with distance 21 

15. Variation of void fraction with distance 21 


JPL TECHNICAL REPORT 32-987 


v 


Contents (contd) 


Figures (contd) 

16. 50-in. nozzle and injector 22 

17. Injection orifice details 23 

18. Experimental nozzle in operation 24 

19. Exit jet at 1 .0-/xs exposure 25 

20. Water injection pressure drop 26 

21. Nitrogen injection pressure drop 26 

22. Freon injection pressure drop 26 

23. Comparison of theoretical and experimental pressure profiles ... 27 

24. Comparison of theoretical and experimental variation of 

pressure with mixture ratio 21 in. upstream of exit 27 

25. Comparison of theoretical and experimental variation of 

pressure with mixture ratio at exit 27 

26. Comparison of theoretical and experimental flow rates at 

1 50-psia nozzle inlet pressure 28 

27. Comparison of theoretical and experimental flow rates at 

190-psia nozzle inlet pressure 29 

28. Comparison of theoretical and experimental flow rates at 

230-psia nozzle inlet pressure 29 

29. Comparison of theoretical and experimental exit velocities 

at 1 50-psia nozzle inlet pressure . . .31 

30. Comparison of theoretical and experimental exit velocities 

at 1 90-psia nozzle inlet pressure . . . 31 

31 . Comparison of theoretical and experimental exit velocities 

at 230-psia nozzle inlet pressure . 31 

32. Comparison of theoretical and experimental flow rates 

using Freon 31 

33. Comparison of theoretical and experimental exit velocities 

using Freon 32 

34. 6-in. nozzle and injector 32 

35. Comparison of theoretical and experimental flow rates 

for 6-in. nozzle 33 

36. Comparison of theoretical and experimental exit velocities 

for 6-in. nozzle 33 

37. Comparison of theoretical and experimental flow rates for 

steam-water nozzle of Ref. 28 33 

38. Comparison of theoretical and experimental exit velocities 

for liquid nitrogen nozzle of Ref. 29 34 


JPL TECHNICAL REPORT 32-987 



Abstract 


Two-phase nozzle flow of liquid drops mixed with a gas is analyzed for both 
the one-component and die two-component cases. Drop volume, drop breakup 
and vaporization, variable fluid properties, and solubility of the gas in the liquid 
are taken into account. Computer programs employing the analysis provide gas 
and liquid velocities and temperatures, as well as drop diameter, as a function of 
distance through a nozzle when initial conditions and fluid properties are specified. 
Wall shear and boundary-layer growth are also computed, and an option provides 
the optimum (maximum exit velocity for given length) nozzle contour. Calculations 
for typical fluids illustrate the effect on nozzle exit velocity of nozzle contour and 
length, initial drop size, mixture ratio, and pressure ratio. Nozzle experiments 
with nitrogen-water mixtures and Freon-water mixtures show that actual velocities 
and flow rates agree with the analysis within 5% . 
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Acceleration of Liquids in Two-Phase Nozzles 


I. Introduction 

High-velocity liquids are a source of power for pumping 
(Refs. 1 and 2) and for magnetohydrodynamic (MHD) 
electrical power generation (Refs. 3 and 4). Liquid veloci- 
ties of 300 to 1000 ft/s are employed in these applications, 
and the velocities are produced by accelerating the liquids 
in two-phase nozzles. 

A two-phase nozzle is a nozzle in which a liquid and a 
gas are mixed at high pressure and low velocity and 
expanded to low pressure and high velocity. The gas 
phase may either be the vapor of the liquid being accel- 
erated, in which case the flow is termed “one-component/ 
or a different chemical species from the liquid, in which 
case the flow is termed '‘two-component/' A one- 
component fluid combination of interest for MHD 
power generation is potassium and potassium vapor; 
one-component combinations that have been investigated 
experimentally are water and steam, and liquid nitrogen 
and nitrogen. The fluid combination of greatest interest 
for MHD power generation is the two-component com- 
bination of liquid lithium and cesium vapor. The two- 
component combination which has been most extensively 
investigated experimentally is water mixed with nitrogen 
or air. 


Typical of the nozzle flow conditions for MHD power 
generation are the following, calculated for a 50-in. -long 
cesium-lithium nozzle designed to be employed in an 
MHD system of 300 kW electric output (Ref. 4). Condi- 
tions at the nozzle inlet are: pressure, 137 psia; tempera- 
ture, 1800°F; velocity, 50 ft/s; void fraction, 0.72; liquid 
flow rate, 186 lb/s (180 lb/s Li and 6 lb/s dissolved Cs); 
gas flow rate (Cs vapor), 14 lb/s; and liquid drop diameter 
large enough (say, 0.05 in.) to result in breakup in the 
nozzle. Conditions at the exit are: pressure, 15 psia; liquid 
temperature, 1790°F; gas temperature, 1610 C F; liquid ve- 
locity, 495 ft/s; gas velocity, 687 ft/s; void fraction, 0.96; 
liquid flow rate, 181 lb/s (180 Ib/s Li and 1 Ib/s dissolved 
Cs); gas flow rate, 19 lb/s; and liquid drop diameter, 
0.019 in. The nozzle is converging-diverging, a require- 
ment resulting from the low sonic velocity in two-phase 
mixtures, with a throat diameter of 4.8 in. and an exit 
diameter of 7.8 in. 

It will be shown in this report that the exit velocities of 
the nozzles of interest for MHD power generation can be 
calculated to an accuracy of 5% given the properties of 
the fluids. This accuracy is possible because the nozzles 
employ spatially uniform mixing of the liquid and gas at 
the nozzle inlet and are long and slender, making the 
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one-dimensional equations for two-phase flow closely 
applicable. The only approximation needed is in deter- 
mining drop diameter, and the method adopted for 
calculating this diameter is adequate for 5% accuracy 
in calculated exit velocity. 

Tangren, Dodge, and Seifert (Ref. 5) derived the equa- 
tions for isentropic flow of an immiscible two-component 
mixture with equal temperatures of the liquid and the 
gas. The first extensive experiments on two-phase nozzles 
were reported by Elliott (Ref. 1), who found that actual 
nozzles achieve 80 to 90% of the isentropic exit velocity. 

The differential equations for two-phase flow with 
velocity and temperature differences between the phases 
were derived by Kliegel (Ref. 6). Those equations were 
successfully applied by Kliegel and others to gas-solid 
flows in rocket motor nozzles, a situation in which the 
particle sizes are known. The fact that particle sizes can 
also be determined for a gas-liquid nozzle was pointed 
out by Crabtree (Ref. 7). He showed that the flow con- 
ditions of the drops in a two-phase nozzle are such that 
the drops will break up early in the acceleration process, 
and he suggested calculating the resulting diameter from 
the Weber number = 6 criterion, which correlates the 
data for atomization by gas streams within a factor of 
about two. An uncertainty in drop diameter of this mag- 
nitude causes only about a 3% uncertainty in nozzle exit 
velocity. The We = 6 criterion was incorporated into a 
two-phase nozzle computer program by Netzer (Ref. 8), 
and the computed exit velocities agreed within 3% with 
measured values in a short air-water nozzle, where 
friction would be expected to be small. A long nozzle, in 
which friction would be expected to be significant, gave 
velocities 25% below the prediction. 

Netzers analysis employed equal phase temperatures, 
and this approximation was eliminated in an analysis by 
Elliott (Ref. 9) which accounted for both the velocity 
difference and temperature difference between the 
phases, while again employing the We — 6 criterion for 
drop breakup. The analysis of Ref. 9 gave nozzle exit 
velocities that agreed within 5% with measurements on 
a large nitrogen-water nozzle. 

This report describes the final version of the compu- 
tation method reported in Ref. 9. The effects added were 
friction, mass transfer between phases through vapori- 
zation and condensation, variable properties, solubility 
of the gas in the liquid, and finite vapor pressure of the 
liquid. The analysis retained the basic assumption of one- 

2 


dimensional flow. A procedure for determining the opti- 
mum nozzle shape was added, and the one-component 
flow case was analyzed. 

II. Analysis 

The problem analyzed is illustrated in Fig. 1. A spa- 
tially uniform one- or two-component mixture of liquid 
drops and gas enters a nozzle at high pressure and low 
velocity and expands to low pressure and high velocity. 
The objective of the analysis is to determine, for a speci- 
fied pressure profile p(x), the drop diameter D and the 
temperatures T g and Ti, velocities V g and V h and flow 
rates ih g and ffp of the gas and liquid phases, respec- 
tively, at each station in the nozzle, given the initial 
values of D, T g> T if V g , V u the total flow rate, and the 
properties of the fluids. 

The five relations employed to compute the five un- 
knowns D, T g , T h V g , and V t are (1) the momentum 
equation for the mixture, (2) the energy equation for the 
mixture, (3) the drop drag equation, (4) the drop heat- 
transfer equation, and (5) the drop breakup criterion. 



A 

Fig. 1. Two-phase nozzle flow geometry 
and nomenclature 
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Solubility and vapor pressure relations provide the flow 
rate ratio rh g /rhi. 

A. Two-Component Flow 

1. Assumptions . The assumptions employed in the two- 
component analysis are as follows: 

1. The liquid is uniformly dispersed as spherical drops, 
all of the same diameter. 

2. The drops break up to limit the Weber number to 6. 

3. There are no external forces acting on the two-phase 
mixture other than pressure and wall shear, and 
there is no heat transfer to or from the mixture. 

4. The flow is one-dimensional. 

5. The drops are large enough for the surface curvature 
to have negligible effect on the vapor pressure of 
the liquid and for the surface energy to be negligible. 

6. The drops are isothermal. 

7. The gas mixture obeys the additive-pressure law. 

8. The partial pressure of the predominantly liquid 
component is given by Raoult’s Law. 

9. The concentration of the predominantly gaseous 
component in the liquid is given by Henry’s Law. 

10. The volume of the liquid solution is equal to the sum 
of the volumes of the pure liquids. 

Assumption 1 restricts the analysis to nozzles having 
spatially uniform injection of the liquid into the gas and 
operating at gas-to-liquid volume ratios greater than unity. 
Assumption 2, the drop breakup criterion, states that drop 
diameter is limited to a value D for which We = p g V\ D/£t 
— 6. Thus, 


D 


max 


12qr 

p*v] 


(1) 


where p g is the gas density, V s is the slip velocity V g — V h 
and a is the liquid surface tension. The form of Eq. (1) is 
physically reasonable in that the Weber number is pro- 
portional to the ratio of stagnation pressure p g V 2 8 / 2 to 
surface tension pressure 4<j/D. Hence, a drop would be 
expected to flatten and break up at a sufficiently high 
value of We. This has been verified experimentally (Refs. 
10 and 11) and the critical Weber number found to be 6, 
within a factor of 2. An additional restriction is that for 
actual breakup to occur, the time spent at a Weber num- 
ber exceeding 6 must be longer than the natural period 
of oscillation of the drop, tt ( piD 3 /o-) 1/5, /4, where p t is the 


density of the liquid (Ref. 11). As shown later, this re- 
quirement is met only in two-phase nozzles longer than 
about 10 in., and Assumption 2 may cause the analysis to 
overestimate the exit velocity by increasing amounts as 
the nozzle length decreases below 10 in. 

Assumption 3 excludes magnetohydrodynamic and 
mechanical body forces. The exclusion of wall heat trans- 
fer is correct for the insulated nozzles of interest for 
power systems. 

Assumption 4 is closely met in practical nozzles since 
good performance requires small wall angles, large throat 
radius of curvature, and uniformly distributed injection 
of the fluids at the nozzle entrance. 

Assumption 5 is valid for the drop sizes of 0.001 to 
0.010 in. produced by the Eq. (1) breakup criterion. 
Assumption 6 is valid because of the rapid internal circu- 
lation in drops (Ref. 12). Assumption 7 introduces negli- 
gible error in most cases of practical interest since the 
vapor pressure of the liquid is small and need only be 
evaluated approximately. 

Assumptions 8, 9, and 10 are either valid, or cause little 
error, for fluids of low miscibility, which are the fluids of 
interest for power systems. 

2. Derivation of equations for ffee-stream flow 

a. Continuity. iTef^ring to Fig. 1, the nozzle flow area A 

equaLto the gas flow area m g /p g V g plus the liquid flow 

area rhi/piV t . Thus, 

A - m g ( — H (2) 

V Ps y 9 Pi y i ) 

where r is the mass mixture ratio rhi/m g . 

b , Momentum. By Assumption 3, the only force acting 
on the free-stream flow is that due to the pressure gradient. 
If M is the momentum flux at flow area A, the change in 
momentum flux across pressure increment dp is 

dtii = —Adp (3) 

The momentum flux can be written as the sum of the 
momentum fluxes of the gas and liquid. Thus, 


i. 1 = m g V g + ThtVj (4) 

If the flow were allowed to continue at constant pres- 
sure, V g and V i would become equal to each other at the 
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mass-weighted mean velocity V. Since, for this process, Substituting Eqs. (11)— (13) into Eq. (9), the differential 
dAI = 0, the value of V is given by momentum equation is 


(rh g + rh ( )V = m,V g + m ; V t 


( 5 ) 


or 


— V 

y r- 


rV, 


( 6 ) 


1 + r 

Thus, the momentum flux can be written 

M = (m g -f rhi)V (7) 

Since m 0 + vi L is constant, the change in momentum flux is 

did = (rhg + m t )dV (8) 

Substituting Eqs. (8) and (2) into Eq. (3), dV can be written 


dV = 


~ 1 + r(p,V, + Pl V>) dp (9) 


The slip ratio is defm»d as 

s = VJV = (V,-V,)/V 


( 10 ) 


This equation can be combined with Eq. (6) to give V g and 
Vi in terms of Vl 


v ' = ( 1 + ttt) ? = “ v 
Vi = (i - ttt) p 


= bV 


( 11 ) 

(12) 


The gas density can be expressed as 
Ps = Wfp/RTg 


2 VdV = dV" = 


+ < 14 » 


A semi-integrated form of this equation can be written 
which permits a larger step size in numerical integration. 
The quantities a and b are slowly varying because s is 
typically only 0.1 to 0.3 and slowly varying. The quantities 
r, T g , W g and p t are also slowly varying. Integrating 
Eq. (14) over a pressure increment A p, for which a , b , r, 
Tg, Wg, and pi are constant to within the desired accuracy, 
the change in V 2 is 




f 2 l 

f RTg 

+ _L_\ 

/ 1 + r\ 

I A p 
f v- — 

v aW g p 

b P i ) 


AV 2 = 


All quantities other than pressure, can be taken outside 
the integral and evaluated at their mean values (denoted 
by subscript m) corresponding to the mid-interval pres- 
sure p. Thus, 


aV 2 = - 


X 


1 + 


RTa 




A P 

2 -L + 


A p 


A p 


P b mP l 


dp 


m I 

p — r 


(16) 


Performing the integrations, 


AV 2 = 


_ 2 / RT *m 

1 T I'm \ &mWg m 


T P + _L 

w a „ 10ge p - Ap/2 b m p lm 


(17) 


(13) Equation (17) is the final form of the momentum equation. 


where W ff is the effective molecular weight of the gas 
mixture and R is the universal gas constant. Equation (13) 
is the definition of the effective molecular weight W g , 
which is the quantity that gives the actual gas density 
when substituted in Eq. (13). 


c. Energy. The enthalpy change of the mixture between 
state 1 (the beginning of pressure interval A p) and state 2 
(the end of the interval) can be evaluated in two steps: 
(1) phase change at p u and (2) change to p >, Tg 2 , 

Ti 2 , at fixed composition. 


4 
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The enthalpy change for step 1 is 

( enthalpy required \ 
to vaporize and \ 
heat unit mass of I 
A from to T ffl / 

/enthalpy required \ 

+ / amount of B\ ^ / to vaporize and \ 
^vaporized J l heat unit mass of I 
\B from 7\ to T Sl j 

/kinetic energy 

/ amount of A and\ ^ / required to 
\B vaporized J l accelerate unit mass 
\from V lr to V 9l 


or 


AHx = (m„ g2 - Tha Si ) [L ttl + c, Vi (T ffl - T h )] 

+ {m bg2 - m brji ) [L», + 

+ (rh g 2 -m Sl )(Vlr V h V 2 < 18 ) 

where L and c are latent heat and specific heat, respec- 
tively. 

Introducing more compact notation, 


ah,- Arh^ (L ai + Ca 9i *iT) 

+ Am,,(I,., + c., i S,r)+-i=!|£- (19) 


The enthalpy change for step 2 is evaluated from the 
temperature, pressure, and velocity changes, with proper- 
ties evaluated at mean T and p for the interval. 

A H 2 = m, 2 [c, m (T S2 - T g J + -j- (V; 2 - V 2 ti )] 


+ ™i 2 [ c i m ( T h~ T h) + 




P2 Pi 

pi m 


( 20 ) 


= AT <> + 

+ m l2 (c lm AT i + -^- + ^ (21) 


By Assumption 3, no work is done by the free-stream 
flow and no heat is transferred to it. Hence, 


A H l + AH 2 = 0 


(22) 


Substituting Eqs. (19) and (21) into Eq. (22) and solving 
for AT* gives the energy equation for the mixture: 


+ 


Am, 


^r + ^r( L - + c --. s ‘ T ) 

+_s r( L '' +c *>' 8 ' T )] (23) 


d . Drag. Although no force other than pressure acts on 
the free-stream flow as a whole, a drag force exists be- 
tween the phases. Hence, a second momentum equation 
must be written using as the control volume the boundary 
between the phases. 

The two forces acting on each liquid drop are the 
buoyancy due to the pressure gradient and the drag due 
to the relative gas velocity. The sum of these is equal to 
the mass times the acceleration of the drop. Thus, for a 
single drop. 


( dynamic pressure of\ / drag \ x ( ^ ronta ^ area \ 
relative gas flow J \ coefficient J \of drop J 

( volumeX / pressure \ __ / mass \ ^ / acceleration^ 
of drop J ^gradient J \^of drop J yof drop J 


or 




7tP 3 dp 
6 dx 


= (24) 


The absolute value sign in the first term makes the 
drag force positive when V g > V; and negative when 
V, < V,. 


Solving Eq. (24) for dV j, 


nr _ | s j sV 2 C D dx dp 

4 Pl ViD p{V : 


(25) 


JPL TECHNICAL REPORT 32-987 


5 



Differentiating Eq. (12), dV i can also be expressed in 
terms of s, r, and V. Thus, 


evaporative cooling due to the latent heat supplied to the 
liquid vaporized. The convective cooling is 


dVi = bdV + V < 26 > 

Solving for ds , 

, _ b (1 + r) dV + sdr _ (1 4- r) dVj ^7) 

~ v 1 -hr v 


-dQ c = hA d N(T t - T g )dt (32) 

where h is the heat-transfer coefficient, A d = ttD 2 is the 
surface area of a drop, and dt = dx/V 1 is the time required 
to traverse dx. Thus, 


Substituting dV 1 from Eq. (25), noting that dV = 
dV 2 /2V y using Eq. (12), and writing for a finite increment, 

__ b m (l 4- r m ) AV 2 (1 4- r m ) A p s m Ar 

2V 2 b mPl V 2 1 + fm 

m l m m 

Zpv m \s m \s m C r>m (l + r m )\x 

4 b mPlm D {2 *> 


This is the drag equation employed when x is specified 
as a function of p. In a later section, the problem will be 
changed to find the optimum value of As and then the 
corresponding value of Ax. Solving Eq. (28) for Ax yields 
the required alternative equation: 


Ax = 


4 D 

$m S m C Dm V’ 


Ap + 


b 2 m p> m *v* 


bm V 2 m / s M Ar 
1 + T m \ 1 + Tm 



(29) 


e. Heat transfer. Although no heat is transferred to the 
mixture as a whole, heat transfer exists between the phases. 
Hence, a second energy equation must be written using 
as the control volume 4h e boundary between the phases. 

The work dW done on the liquid is that due to drag 
by the gas. (Only work done by shear or shaft forces is 
included in dW when writing the First Law for a control 
volume, as shown in Ref. 13, p. 39.) Multiplying Eq. (24) 
by the number flow rate of drops N — 6rhi/7rD 3 p ly the 
drag force F d on that quantity of liquid is 


F d = A p, | V a | V s C d wD 2 = ^ + m t v , dVl 


dx 


The work done on the liquid is 


— dW = F d dx — m 


. (dp dVj \ 


(30) 


(31) 


The heat dQ transferred from the liquid is made up of 
two parts: (1) the convective cooling due to the tempera- 
ture difference between the liquid and gas and (2) the 


-dQr = 


6hrhi(Tg — Ti)dx 
DpiV 1 


The evaporative cooling is 

— dQ v = L a drh ag 4- L b dm bj 


(33) 


(34) 


The change in enthalpy of the liquid over the pressure 
increment dp is 

dH = m [c,dTi + & + (35) 

Substituting Eqs. (31), (33), (34), and (35) into the 
steady-flow energy equation dQ — dW — dH, the result is 

6 ^Dp!vf X ~ ~~ ~ L bdm bg = miCidTi (36) 

where ST = T g — Ti. 


Writing for a finite interval, the final form of the drop 
heat-transfer equation is 


ATi 


i r 6fts m r,Ax _ La 

C *m _ D P l m V ^m * mi m 


J 


(37) 


Equations (1), (17), (23), (28), and (37) are the five equa- 
tions that must be solved simultaneously to obtain the 
values of the five dependent variables D, T g , T h V g , and 
Ve as a function of the independent variable p. To carry 
out the solution, all quantities in the equations must be 
expressed in terms of these six variables. The additional 
equations required will be developed next. 

f. Phase properties. By Assumption 9, the mole fraction 
of component A in the liquid is proportional to the partial 
pressure p a of component A (Henry’s Law, Ref. 14, p. 
181). Thus, 
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m ai fW H 


Th ai /W ai + ( rhi ~~ iha l )/Wi >l 


= Hp a (38) 


Substituting 6, from Eq. (45) then gives expressions for 
m ag and m bg in terms of molecular weights and pressures: 


where m ai is the flow rate of dissolved component A, W aj 
and W),, are the molecular weights of liquid A and B, 
respectively, and H is the mole fraction dissolved per unit 
pressure of A. 


m ag = m,j 


. w agPa 

W sP 


. _ . W K p b 

m tg Wsp 


(47) 

(48) 


Solving for = « yields 

(W ai /W tl ) Hp a 


a 


I + (W H /W bl - 1) Hp„ 


(39) 


From Eq. (48), the value of P — th bg /ihg, the fraction of B 
vapor in the gas, is 


By Assumption 8, the partial pressure p b of component 
B is equal to the vapor pressure p bg of pure B multiplied 
by the mole fraction of B in the liquid (Raoult’s Law, 
Ref. 14, p. 101): 

Pi = (1 - ffp«K (40) 

But by Assumption 7, 

p = Pa + p b (41) 

Combining Eqs. (40) and (41) and solving for p a yields 

V ~ P» „ 


1 w g p 


(49) 


The mass ratio of component B to component A, 
r c = mb/rha, can be written 

= + mi, _ m t - rh a; + m ^ 

Tc m ai + m ag m g - m bg + m H K 

Solving Eq. (50) for r = thi/m g and employing the 
definitions of « and p, the mass ratio of liquid to gas is 


r = 


fc ~ (1 + rjp 
1 — (1 + r c )a 


(51) 


Pa 


1 - HPb o 


(42) 


By Assumption 10, the volume flow rate of each com- 
ponent in the liquid is 


By Assumption 7, the mass flow rate of each component 
in the gas phase can be written 


c«, = 


m a , 
Pa l 


(52) 


”S = 

Wa gPa 

(«) and 



RTg 



rh »g = 

w bgPb 

RTg Vg 

(44) 

m bl mi — m ai 

Vbj = — 

P*>i P b i 

(53) 


where W„ g and W bg are the effective molecular weights 
of gaseous A and B, respectively, and v g is the gas volume 
flow rate. 


The gas flow rate from Eq. (13) is 


where p H and p b{ are the densities of liquid A and B, 
respectively. 

From Eqs. (52) and (53) and the definition of a, the 
liquid density is 


• _ w gP 

Wig RT ^ 


(45) 


pi = 


(a/p n ) + (1 — a)/p bl 


(54) 


Adding m ttg and m bg and equating to m B yields 
W 0 p a + W b p b 


The specific heats of the phases are the mass-weighted 
averages (Ref. 15, p. 494). 




(46) 


Cg (1 P) C ag + pC bg 


(55) 
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Ci = a c ai + (1 — a)c bl (56) 

The liquid viscosity is taken as the mass-weighted average: 
fit = a {L ai + (1 — a) ji bl (57) 


g, Coefficients. The Reynolds number of the slip flow 
between the gas and a drop is 


Re = 


p g DV [ s 

fig 


(64) 


The relation adopted for the viscosity of the gas mix- 
ture is Wilke's equation (Ref. 15, p. 494): 

_ ^9 , ^ /tO\ 

^ 1 + (X t /X a ) <f > «» 1 + (X a /x b ) <t> ba (00) 

where X b /X a is the mole ratio given by 


The relationships employed for computing the drop 
drag coefficient C D from Re are plotted in Fig. 2. For 
Reynolds numbers less than 0.1, the drag coefficient is 
given by Stokes Law: 

24 

C D = — Re < 0.1 (65) 


_ ( i \ 
x a 

and <f>ij is given by 

. _ [i + jitMT (w if /w it )^ 
2 V2 [1 + W i ,/W„]'* 


(59) 


(60) 


The relation adopted for the thermal conductivity of 
the gas mixture is the Lindsay and Bromley relation 
(Ref. 15, p. 494): 


For 0.1 < Re < 2 X 10 J , the Co data for solid spheres 
are represented by Stonecypher’s least-squares fit (Ref. 16, 
p. 3) to Perry’s tabulation (Ref. 17, p. 5-61): 

log, C D = 3.271 - 0.8893 (log. Re) + 0.03417 (log,. Re) 2 

+ 0.001443 (log,. Re) 3 0.1 < Re < 2 X 10* 

( 66 ) 


For Reynolds numbers greater than 2 X 10 C„ is 
taken as constant at 0.4569, the Eq. (66) value for Re = 
2 X 10k 


k a = 


+ 


1 + (Xt/Xa) B„ 1 + (Xa/Xt) B 


(61) 


where 


B a = 


1 

4 


{[ 


1 + 


(Wt v*/l + S a /T s YV*\ 2 
WaJ \1 + S b /T g )] / 


i + VsTsl/T g 

1 + S a /T s 


(62) 


and 


B b = 



W a \*/l + s 6 /TAp 
W bg ) \1 + SJTjj / 


i + VsTs^/r, 

1 + S b /T s 


(63) 


The quantities S a and S 6 are the Sutherland constants for 
the components, and may be roughly approximated by 
1.5 T Bi where T B is the normal boiling temperature. 


The Co values for solid spheres are only a rough 
approximation to the values that can be encountered with 
liquid drops, because the drops can deform into other 
shapes. The data bands for three sets of drop experiments 
are included in Fig. 2. Ingebo (Ref. 18) measured drag 
coefficients for drops of 0.001- to 0,005-in. diameter rapidly 
accelerating in air streams; at the highest Reynolds num- 
bers of the experiments, the C {} values were only a third 
of the solid-sphere values. Hughes and Gilliland (Ref, 19) 
surveyed the information on drops falling in air and found 
that the drag coefficient became higher than the solid 
sphere values at Reynolds numbers above about 200. 
Rabin, Schallenmuller, and Lawhead (Ref. 20) observed 
the acceleration of fuel drops in a shock tube and found 
Co values as high as 2.0 at a Reynolds number of 10 l , 
where the solid sphere value is 0.5. The Reynolds num- 
bers for the two-phase nozzles discussed in this report 
were all in the range of 100 to 4000, with the higher 
values at the inlet and the lower values at the exit of the 
nozzles. 

Thus, there are two factors which are expected to make 
the drag on the liquid phase larger in practice than 
assumed in this analysis. First, the drop size in practice 
will have the Weber number = 6 value as the upper 
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limit, whereas the analysis employs that value for the 
mean size. Second, the drag coefficient in practice, at least 
at the higher Reynolds numbers, is larger than the solid- 
sphere coefficient assumed in the analysis. 


From Ref. 21 (pp. 414-415), 
coefficient is taken as 

the drop heat-transfer 

s* 

II 


Re < 1.0 

(67) 

b c y p v V | s | ^ 

'2.2 0.48 > 

. Re \/Rey 

) 1.0 < Re <25 

(68) 

, 0.37 k. Re 0 

h ~ D 

.6 

Re > 25 

(69) 


These equations complete the analysis of frictionless 
two-component flow with a prescribed pressure profile 
p(x). 

3. Optimum pressure profile. The optimum pressure 
profile is the one that maximizes the mean exit velocity 
(which is also the velocity attained by the liquid after 
exhausting into a constant-pressure region) for a given 
nozzle length L. That is, it is desired that V e , or, for greater 
mathematical convenience, VJ, meet the requirement 


v - =v ° + Lw+- 

subject to the constraint that 

fpe dr 

x e = J ~fl^dp — fixed 


maximum 


(70) 


(71) 


Since the equations of the preceding section can be 
solved by specifying the slip s directly as the independent 
variable instead of calculating s from the pressure profile 
through Eq. (28), the integrand of Eq. (70) can be con- 
sidered a function of s and p ; namely, from Eq. (14), 

dV 2 2/1 r \ 

W ~ - TT7(-^7 + TpJ = F ( s > v) (72) 

From Eq. (29), the integrand of Eq. (71) is 


dx 


4D 


dp 3/3, |sjsC D V 2 [ * ' '2 dp 

+ ML(-S_d L _ds\- I 
1 + r \1 + r dp dp) J 


1 + 


b 2 pi dV 2 


Numerical estimates show that dr/dp and ds/dp are 
small enough in a two-phase nozzle so that the last term 
in Eq. (73) can be neglected. Hence, dx/dp is approx- 
imately 


dx 


4D 


-z— ( i + = G(sn) 

dp 3 P! ,\s\sC d V 2 \ 2 ) 

(74) 


Thus, the problem is to find the function s(p) for which 

e 

F (s, p) dp ~ maximum (75) 

subject to the constraint that 


/’ 


/' 


G ( s , p) dp = fixed 


(76) 


This is the isoperimetric problem of the calculus of 
variations. The desired maximizing function, designated 
s 0 {p), is found (Ref. 22, pp. 216-218) by forming the 
function 


F* = F- \G 


and solving 


dF* 

ds 


= 0 


(77) 


(78) 


where A is a constant, the Lagrangian multiplier, which 
has the units of acceleration here. It is a constant which 
specifies the length of the nozzle for given fluid properties 
and inlet conditions; the method of selecting it will be 
discussed in Section IV. 

Substituting Eqs. (74) and (77) into Eq. (78), the opti- 
mum slip is the solution of 


0F 4Da 9 


[Rk(-¥)] = « 


3* 3p, V* ds 

Combining the terms in dF/ds, this becomes 

A _ 2Dxb 2 P i \ 0F 2DA piF 0b 2 

V 3 p„7 2 \s\sC D ) ds 3 p,V 2 I s j s C D ds 


(79) 


- JSL ( i + M\i/ M = 0 

3p,V 2 \ 2 /0s^|s|sC O/ / 

(80) 
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Evaluating F from Eq. (72), performing the indicated 
differentiations, and solving for s, the relation for the 
optimum slip (for frictionless flow) is 


Ape’* ( fgs 0 + eC> D y 
| Ap I \ 2g (d — 1) / 


(81) 


where 


and 


c — 

d = 

e = 

f = 

g = 


4D\ h 2 pi 
3 Pff V 2 Cn 

fc 2 pi 

a 2 p g 


h 2 pi , br 

ap g (1 + r) 1 + r 



r 

(1 + r) 2 


-1 


(82) 

(83) 

(84) 

(85) 

(86) 


is nearly isothermal, so that p„ is proportional to the 
pressure p, then Eq. (81) reduces to 

= const X (-L) (90) 

To utilize Eq. (90), the V(p) relationship for isentropic 
flow can be substituted to give the optimum slip s 0 as a 
function of pressure for a given initial slip. The optimum 
slip can then be substituted in Eq. (74) to give the opti- 
mum dx/dp from which x(p) can be calculated. The 
nozzle length is determined by the initial slip assumed. 

4. Wall shear and boundary layer . The momentum 
thickness 0 of a boundary layer is the thickness of free- 
stream flow which has a momentum flux equal to that 
by which the momentum flux of frictionless nozzle flow 
exceeds the momentum flux of the real flow for the same 
mass flux (Ref. 24, p. 5). For a two-phase nozzle, the 
momentum flux of the frictionless nozzle flow is that 
given by Eq. (7): 


M = rh t V 


(91) 


C' D = 1.0 


He < 0.1611 


(87) 


The mean mixture density corresponding to the mean 
velocity V is 


C' D = 1.1107 + 0.06854 log. Re 
+ 0.004329 (log. Re) 2 

0.1611 < Re < 4.709 X 10 3 

( 88 ) 


P' = 


rht_ 

AV 


Pi 

1 + r a 


(92) 


where r a is the ratio of gas flow area to liquid flow area 
PlVl/rp g V g . 


C' D = 2.0 Re > 4.709 X 10 3 (89) 

A similar analysis was employed by Marble (Ref. 23) 
to find the optimum nozzle shape with Stokes flow. 

In the computer programs the (d — 1) term in Eq. (81) 
is replaced by d , with negligible error at the small slip 
values of a high-performance nozzle; the substitution 
avoids numerical difficulties in calculating s 0 for the initial 
pressure step in some nozzles, where the slip could 
initially be high enough to give d = 1, and the results 
should be considered invalid in any part of a nozzle 
where d is less than about 5. 

Equation (81) can be simplified to give a rule for find- 
ing the optimum nozzle contour for many cases of prac- 
tical interest. If r and pi/p g are large, s 0 is small, and 
Re does not vary greatly, the quantity in parentheses in 
Eq. (81) is nearly constant. If D is constant and the flow 


From the definition of the momentum thickness, the 
value of 0 at a station where the nozzle wall radius is y 0 
is given by 

M - M, = 2 wij„e = 2t tij„p’V-6 (93) 

where M / is the momentum flux of the real flow with 
friction. 

Equation (93) is identical to the definition of 8 for 
single-phase flow (Ref. 24, Eq. 11), with the single phasep 
replaced by p' and the single-phase V replaced by V. 
The skin-friction coefficient can be defined using the same 
quantities. 
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where t w is the wall shear. It will be shown that a valid 
Cf value can be provided. 

The boundary-layer momentum equation (Ref. 24, 
Eq. 22) then becomes 



For a W-power profile, the velocity thickness, from 
Eq. (12) of Ref. 24, is 

s = - j-e (ioo) 

A convenient relation for C/ z as a function of Rd is the 
Shultz-Grunow relation (Ref. 15, p. 147) which, with the 
aid of Eq. (6-35) of Ref. 15, can be written 


(95) 

where 8* is the displacement thickness, i.e., the distance 
the wall must be moved outward to give the same flow 
rate as with frictionless flow. 


0.208 

(log,, Ra + 0.425) 2 - 584 


( 101 ) 


Comparison of Eqs. (94) and (98) shows that C f can 
be written 


C/ 


rb 

1 + r 


Cft 


(102) 


Assuming a l 4-power velocity profile and no density 
variation (p f ~ const and S' = S in Ref, 24), the shape 
factor 8-/0 is obtained from Eqs. (9) and (12) of Ref. 24. 


/.‘n-fartn* 9 

9 f 6 (y/&y n [i - (y/s)' 7 ’] dy 7 

./ 0 


(96) 


where S is the velocity thickness of the boundary layer. 

Noting that dV can be written dV 2 / 2V, and that 
d(p'V) = d(m t /A) from Eq. (92), the finite-difference 
form of Eq. (95) is 


C f / 8 _ 1 1 \ 

49 = ~f “ - "■ (wr 4VS “ X 44 + ^ A,J • ) 

m 

Wall shear in homogeneous two-phase flow has been 
found (Ref. 25, p. 181) to be equal to that which would 
exist with pure liquid at equal velocity and boundary- 
layer thickness, multiplied by the wetted wall fraction: 


Thus, the final form of the boundary-layer momentum 
equation is 


A0 - 


r m b n 


1 + r m 


Cr,AX — e„. 


• 8A V 2 
7 Vi 


AA 


Ai/A 
y°m ) 


(103) 


Let Va be the mean velocity of the flow including the 
boundary layer. Then, from Eq. (93), 

m ( Va = Mf— ih t V — 2-n-f/,, p'V 2 6 (104) 


Hence, employing Eq. (92), the mean exit velocity 
including the boundary layer is 


Va — V 




(105) 


By the definition of the displacement thickness, the 
flow rate is reduced by the throat displacement thickness 
S? to 



(106) 


Cq v „ Ai _ C fl p, V 2 
2 Pl 1 A ~ 2(1 + r„) 


(98) 


where Cf x is the skin friction coefficient for liquid at a 
Reynolds number of 


= El Yl 8 


(99) 


5. Auxiliary quantities . The foregoing equations pro- 
vide a solution for the basic dependent variables of the 
two-component nozzle flow problem. In this section, ex- 
pressions will be presented for additional parameters that 
are useful in nozzle design and cycle studies. 

a. Equilibrium exhaust conditions. If the nozzle ex- 
hausts into a constant pressure region such as the atmo- 
sphere or the interior of a separator or condenser, then 
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the velocities and temperatures of the two phases can 
reach equality at V and T, respectively. The velocity V is 
that given by Eq. (6). The equilibrium temperature T 
is found by applying Eq. (23) to the constant-pressure 
process: 




i r 

"av; / — 

AVy\ 

A T, = T 

~T, 


£ + r \ C >m* T ‘ + 

t) 


Affjg, 

/«V*\ 

Aril, 


+ 

rh<j 

(— ) + 

ih 9 (L a + c aff ST) 



The outboard half-area is given by 


a _ ^ R%) 

2 COS <f> 


( 112 ) 


The radius from the axis to the flow centerline is 


R<l> = F4 0 — x sin <t> (113) 


Solving Eqs. (Ill) and (112) for R t and R 0 yields 


+ ^r( L * + ^ 8r )] (107) 

where unsubscripted quantities are those at the exit, 
barred quantities are those after equilibrium is reached, 
barred A’s are changes between the two, and sub ms 
designate quantities evaluated at the mean temperature 
between the two. 


Ri = 4^^ ( 114 ) 

«• = (115) 


Solving for A Tj, 


A Ti 


1 " 

mmr C Q 

’ Ci " 


— A V- + r AV? 

aT, + ±- 


+ 


Arft a 


m 


± (L a + c« g ST) 


Am^ 

m a 


(Lb + Cb g 8T) 



(108) 


Setting T — T and V g = V t = V in the phase property 
and other relations yields the remaining equilibrium ex- 
haust conditions. 


b . Nozzle geometry. For a circular nozzle (Fig. 1) the 
radius y 0 is given by 


yo = 4^ ( ioq ) 

The wall angle w 0 is given by 

tan 0,0 = l±x~ ( no ) 

The geometry for an annular nozzle is shown in Fig. 3. 
The inboard half-area is given by 


A _ 

2 cos (j> 


(in) 



Fig. 3. Annular nozzle geometry 
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The wall coordinates are 


and 


Vi = 


R$ — Ri 
cos <f> 


Reynolds number of the liquid film at the mid-area of 
the cone will be 


(116) 


Re f 


_ PiVtf 

pi 


(120) 


Vo 


R 0 — Rtf, 
cos <j> 


The wall angles are 


(117) where V is the film velocity and t f is the collected film 
thickness at circumference (2 ttA) V 2 . For zero void fraction 
of the film, V is given by 


tan <oi 



(118) 


- PlVtfa rA)* (121) 


and 

tan Wo = 4^ (119) 

AX y 

The annular geometry can be employed in studying 
rectangular nozzles by selecting a large value of and 

a flow rate that gives the desired nozzle height y 0 + JR. 

c. Impingement parameters. For convenience in deter- 
mining the velocity reduction of the liquid due to im- 
pingement on the wall of a separator or condenser, two 
quantities required in the separator analysis of Ref 26 
are computed. 

Assuming that impingement takes place on a cone of 
base area equal to the equilibrium jet area A, the 


Table 1. Conversion, factors employed in 
computer programs 1 


Equation 

number 

Multiply 

By 

1 

Right side 

0.02645 

2, 133 

Right side 

144 

17, 134 

Right side 

4633.04 

23, 108, 130 f 

j V 2 terms 

25,036.47 1 

132, 136, 137 ) 

1 Ap 

5.4039' 1 

28, 29 

Ap 

4633.04 

37 

h m 

3600' 1 

64 

Right side 

300 

67, 69 

Right side 

12 

68 

Right side 

3600 

82 

Right side 

1 2 -1 

98 

Right side 

4633.04 1 

122 

Right side 

43,200 

123 

Right side 

1 44 _1 

* Based on 32.1739 lbm-ft/lbf-s 2 , 778.161 ft-lbf/Btu, 0.03281 ft/cm, 
2.248 X 10-° Ibf/dyn. 


Hence, the mean film Reynolds number is 

Re, = — ’K— (122) 

U 8*A)“ 

The liquid velocity reduction depends on the ratio of 
the mass flux for all-liquid flow to the actual mass flux. 
This ratio is 



which is 2/C/ times the Q defined in Ref. 26. 

6 . Units . The units employed for each quantity are 
listed in the Nomenclature. Computer input and output 
quantities are in those units. Those equations requiring 
conversion factors, and the values of the factors, are 
listed in Table 1. 


B. One-Component Flow 

1. Assumptions . Assumptions 1-6 of the two-component 
analysis are retained. The additional assumption is made 
that both phases are at the saturation temperature corre- 
sponding to the local pressure. This is realistic in that 
evaporation of the liquid and condensation of the vapor 
strongly drive the temperatures of both phases toward 
the saturation value. 

2. Derivation of equations for free-stream flow 

a. Continuity. Equation (2) is still valid for this case. 
One notation change is made, however. The mixture 
quality X g — m g /m t is employed in place of the mixture 
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ratio r = rhi/rhg throughout the one-component program, d . Drag. Equations (28) and (29) correctly express the 

and r is everywhere replaced by drag relation for this case. 


r = 


L. L 

A, 


and changes in r by 


A r ~ — 



(124) 


(125) 


b. Momentum. Equation (P; is valid for this case 
except that, in accordance with the equal-temperature 
assumption, T &m is simply the mid-interval saturation 
temperature T m . 


c. Energy. The total enthalpy of the vapor phase is 

H g = m, (h, + (126) 

where h g is the saturated vapor enthalpy, and the total 
enthalpy of the liquid phase is 


e. Heat transfer. Heat-transfer considerations are elim- 
inated by the assumption of equal temperatures for the 
two phases. The one-component case reduces to four 
equations (1, 17, 28, and 130) in the four unknowns, 
D, V„ V,, and X,. 

3. Optimum pressure profile. Equation (81) is valid 
for the one-component, as well as the two-component, 
case. 

4. Wall shear and boundary layer. Again, the equations 
are the same as for the two-component case. 

5. Auxiliary quantities. The energy equation for attain- 
ment of equilibrium at constant pressure is 

dig (h g H — ^ + rhi(lii + = m t V 2 + mji g + mlhi 

(131) 


H l = mi(c l T + ^-+x) (127) 

The energy equation for the mixture is 

a[ m, (h, + + rhi ( c{T + + -^-)] = 0 

(128) 


or 


Am, (h Sm + -^) + rhg m (a h, + -^) 

+ Am< (c m T + ( 129 ) 

+ rhi m ( c lm AT + ^L+^l)=0 

Dividing by rh t and rearranging gives the quality change 
across the pressure increment. 


AX, = - 


V 2 — V 2 
L m + 11m lia 


x*. (a h t + 


+ 


(i - *.„) ( 


A p A V‘j 

Ar + 'pu + 


¥)] 


(130) 


Dividing by m f , setting hi = hi, and rearranging gives 
_ / v- \ V- V- 

x ' = x '(i + if) + (1 - x ->zi-w 

(132) 

6. Units. The conversion factors for the one-component 
equations are listed with the two-component factors in 
Table 1. 


C. Isentropic Flow 

In the limit as drop size approaches zero, the phase 
temperatures and velocities approach equality and the 
nozzle flow becomes a reversible, isentropic process. The 
corresponding exit velocity is the maximum attainable for 
any given inlet conditions. 

1. T wo-component flow 

a. Continuity. With V g — V h Eq. (2) reduces to 

A = J rr(—+ —) (133) 

V \p s pi) 

b. Momentum. With s = 0, Eq. (17) reduces to 


AV- 


2 / RT m . p + Ap/2 r„,Ap\ 

1 + P-AP/2 Pi m ) 

(134) 
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c. Energy. With V g = Vi and T g — Ti, Eq. (23) reduces 
to 


i r aV 2 

AT = -— -^ + r 2 [c l AT + 


+ 


Am- 


L 2 


( e 


A p 

P'm 


+ 


AV 2 ' 


A m b/i 


m 


V-2 


La, + L b , 


m Q 


(135) 


Collecting terms in A T, this becomes 
1 


AT = 


+ r > c ln 


(1 + *) ^ 


Ph, 


+ 


L„, im„ ( + L bl \m bq 


m 


(136) 


@2 


2. One-component flow 

a. Continuity. Equation (133) applies to this case. 

b. Momentum. Equation (134) applies. 

c. Energy. With V g = V t , Eq. (130) reduces to 


A Xg — 

+ (i - *..) ( 


+ 


aV- 


.rpj_ A p AW 

c i nt AT d 1 5” 

P'» 2 


(137) 


III. Computer Programs 

Computer programs were written for (1) isentropic 
two-component flow, (2) real two-component flow, (3) 
isentropic one-component flow, and (4) real one- 
component flow. (“Real flow” refers to flow with slip.) 
The real flow programs compute both the frictionless 
free-stream flow and the friction loss due to the boundary 
layer. The prescribed quantities are the pressure, drop 
diameter, and phase velocities and temperatures at the 
nozzle entrance; the exit pressure; the number of pressure 
steps desired (for balancing accuracy against computing 
time); and the fluid properties as a function of tempera- 
ture or pressure, or both. 


The programs begin by storing the tables of fluid prop- 
erties and computing the nozzle inlet conditions. The 
programs then proceed half a pressure step at a time. 
At the middle of each pressure interval, the changes in 
quantities across the interval are computed using the 
properties interpolated from the tables for that pressure, 
and for the existing temperature in the case of those prop- 
erties tabulated as functions of both p and T. The change 
in slip is found if p(x) is specified, or a new optimum slip 
is found if the nozzle is being optimized. At the end of 
each pressure step, the flow conditions are updated and 
initial conditions are determined for the next step. The 
drop size is reduced to D max (Eq. 1) at that time if the 
Weber number exceeds 6 (D exceeds D max ) and the 
breakup option has been specified. Flow conditions are 
printed if the pressure is one selected for output. The 
computations continue until the last pressure step has 
been completed, and the flow conditions at the smallest 
flow area encountered are printed as the throat conditions. 

IV. Two-Phase Nozzle Design 

The method of utilizing the computer programs to 
design a two-phase nozzle will now be presented. The 
design illustrated will be that of the 50-in. nitrogen-water 
nozzle employed in the experiments described in Sec- 
tion V. The performance results calculated here will also 
illustrate the characteristics of two-phase nozzles and 
provide values for comparison with the experimental 
results. 

The flow conditions for which the nozzle is to be 
designed are: inlet pressure p 0 = 150 psia (typical of inlet 
pressures for liquid-metal MHD systems); water/nitrogen 
mass ratio r c ~ 40 (typical of MHD systems); inlet tem- 
perature Tg 0 = Ti 0 = 520°R (mean water temperature in 
test facility); inlet velocity V/ 0 = 7.3 ft/s (velocity in 
the 14-in. -diam inlet of the experimental nozzle); total 
mass flow rate rh t — 160 Ibm/s (flow rate calculated for 
the throat area of the experimental nozzle); and nozzle 
exit pressure p e — 14.1 psia (atmospheric pressure at the 
test facility). 


Two options are available for the real-flow programs; 
in the MOP = 0 option, the variation of pressure with 
distance is specified; in the MOP — 1 option, the pressure 
variation is optimized. Nozzle geometry and boundary- 
layer growth can be computed for either a circular nozzle 
(MGEO = 0) or an annular nozzle (MGEO = 1). In addi- 
tion, the liquid drops can be assumed to have either a 
constant diameter (MBU — 0), or to break up to limit 
the Weber number to 6 (MBU = 1). 


The first step is to determine the optimum nozzle 
contour using the program option MOP = 1. This option 
proved useful only for the constant drop diameter 
option MBU = 0. The optimum pressure gradient at the 
inlet is so steep that immediate breakup to fine drops is 
called for by the Weber number = 6 criter ion , while the 
exposure time is too short for that criterion to be valid. 
Suppression of the breakup until the exposure time 
reaches the natural period of oscillation of the drops 
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would have to be added to the programs to permit useful 
results from simultaneous utilization of the optimization 
and breakup options, but this would have no practical 
value because the initial part of the optimum nozzle has 
an unpractically steep wall angle which must be re- 
placed by a more gradual inlet in practice. The optimiza- 
tion option is, therefore, carried out for constant drop 
diameters which bracket the diameter range encountered 
in the final design. The range is 0.01 to 0.05 in. (250 to 
1250 jx m) for this nozzle. 

Computer runs are first made for a series of values of 
the Lagrangian multiplier A covering the range of nozzle 
lengths of possible interest. Figure 4 shows the variation 
of nozzle length L (value of * at p = 14.1 psia) with A for 
the two extreme drop diameters. This figure identifies the 
range A — 10--10* as encompassing the lengths of possible 
interest, i.e., 10 to 1000 in. For A values bracketing the 
range of interest, say 100 and 3000, the variation of nozzle 
radius [/„, normalized by the throat radius is then 
plotted as a function of fractional distance x/L through 
the nozzle, as shown in Fig. 5. The purpose of this plot is 
to determine whether there is any one contour that satis- 
factorily fits the optimum contours for the range of lengths 
and drop diameters being considered. This is seen to be 
possible, in this instance, where the normalized optimum 
contours for two widely differing lengths at both of the 
two drop diameters are closely similar. Any one of 
the contours can, therefore, be adopted as the basis of a 
practical nozzle contour. 

A preliminary guess must now be made as to the opti- 
mum nozzle length to permit modifying the optimum 
profile into a practical nozzle contour. Too short a nozzle 
will have a low exit velocity because of large slip, and 
too long a nozzle will have a low velocity because of 
large friction. Making the guess that the optimum length 
corresponds to A = 3600, the contour for that case, with 
a 0.01-in. drop diameter, is then plotted as shown by the 
dashed curve in Fig. 6. This curve has a nearly vertical 
slope at the inlet. A conical entrance of 20-deg slope is 
probably the steepest that is practical, and the actual 
nozzle inlet is, therefore, made a 20-deg line tangent to 
the optimum contour. The inlet radius, dictated by the 
injector diameter, was 7 in. for the experimental nozzle. 
A conical exit is also chosen to replace the flaring exit of 
the optimum profile, and an angle of 2.5 deg was adopted 
for the experimental nozzle. The resulting practical 
nozzle contour is given by the solid lines in Fig. 6. 
The length of this nozzle is 51.4 in., of which the first 
12 in. is only a transition section between the injector 
and the station where significant acceleration begins. 



LAGRANGIAN MULTIPLIER A, ft/s 2 

Fig. 4. Correspondence between nozzle length 
and Lagrangian multiplier 

The remainder of the design analysis utilizes the pre- 
scribed pressure-versus-distance option MOP = 0, with 
the p(x) profile corresponding to the adopted nozzle 
contour. From the A = 3600 results, the values of pressure 
as a function of distance are plotted utilizing the printed 
x-values for each pressure in the portion of the optimum 
nozzle that was retained. For pressures that fall in the 
conical extensions, the nozzle radius values are read from 
the computer output, and the corresponding x-values are 
calculated for the conical sections. The resulting function 
p(x) is plotted in Fig. 7. 

The next step is to read off a set of values of p and x 
to make up thep(x) input table required in the prescribed- 
pressure option MOP — 0. The values chosen are shown 
superimposed on the curve in Fig. 7. 

Next, the initial drop diameter for the drop breakup 
option MBU = 1 must be selected. Figure 8 shows the 
variation of exit velocity with initial drop diameter for 
three different nozzle lengths. The 50-in. nozzle employs 
the p(x) values from Fig. 7; the other two lengths are 
obtained by multiplying those x values by constant 
factors. The computations show that, with breakup, the 
exit velocity decreases with increasing initial drop diam- 
eter until a diameter is reached above which the velocity 
is constant. This occurs because any larger drop breaks up 
to that diameter. If breakup is prevented (MBU = 0), the 
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NOZZLE RADIUS y 0 ,tn. RADIUS RATIO Y 0 f Yo f 



FRACTIONAL DISTANCE FROM NOZZLE INLET x/L 

Fig. 5. Normalized profiles of optimum nozzles 

velocity decreases indefinitely with increasing initial drop 
diameter. 

The critical initial drop diameter, above which the exit 
velocity is constant, is only 0.020 in. even for a 250-in. 
nozzle. A practical injector with a pressure drop of only 
5-10 psi such as that allowed in an MHD conversion 
system provides few drops this small, and the operating 
regime of interest is, therefore, on the constant-velocity 
portion of the curves. To be sure of being in this region, 
an initial drop diameter of 0.05 in. is used in the sub- 
sequent computations. 

Conformity with the exposure-time requirement of the 
Weber number — 6 criterion is checked at this point. 
Figure 9 presents the variation of drop diameter with 
distance in the 50-in. nozzle. The 0.05-in. drop starts to 



Fig. 6. Modification of optimum shape to 
obtain experimental nozzle 
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Fig. 7. Theoretical pressure profile for experimental 
nozzle and tabulated values employed in 
real-flow analysis 


break up 10 in. from the inlet and completes its breakup 
to 0.011-in. diam at 18 in. from the inlet, 4 in. upstream of 
the throat. The elapsed time from injection can be obtained 
by integrating dx/V i using the computer output values of 
x and V*. The drop diameter as a function of time from 
injection is shown in Fig. 10. The drop diameter is reduced 
from 0.05 to 0.025 in. in 5 ms, and to 0.011 in. in 15 ms. 

The natural oscillation period of the 0.05-in. drop is ^ 

4 ms, and of the 0.011-in. drop 0.4 ms. Thus, the exposure £ 

time is adequate for the Weber number = 6 criterion to ^ 

be valid in the 50-in. nozzle, but probably inadequate in § 

a 10-in. nozzle. ^ 

o 

Finally, the optimum nozzle length can be determined. 
Figure 11 presents the exit velocity with friction included, 

Vs, as a function of nozzle length. For a flow rate of 
160 lbm/s, the exit velocity increases rapidly with in- 
creasing length (because of decreasing slip) until a length 
of about 25 in. is reached. The velocity attains its maxi- 
mum at a length of about 50 in. and then slowly decreases 
because of increasing friction. For 10 lbm/s, the optimum 
length is 12 in. and for 1000 lbm/s, it is about 80 in., 



Fig. 8. Effect of length and initial drop diameter on 
exit velocity, with and without breakup 



Fig. 9. Variation of drop diameter with 
distance from inlet 
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EXIT VELOCITY kg, ft/s 



Fig. 10. Variation of drop diameter with travel time 



Fig. 1 1 . Effect of nozzle length and flow rate 
on exit velocity with friction 


although the cone angles would be excessive in the latter 
case, requiring lengthening the nozzle to perhaps 120 in. 
For infinite flow rate, the friction loss vanishes and the 
velocity increases indefinitely with length; the curve for 
this case is the free-stream velocity V. 



0 10 20 30 4 0 50 60 

DISTANCE jr, in. 


Fig. 12. Optimum and off-optimum pressure profiles 
and corresponding exit velocities 

To illustrate the sensitivity of nozzle performance to 
departure from the optimum contour, Fig. 12 shows two 
pressure profiles which are markedly different from the 
optimum one. The exit velocity V for each profile is 
indicated. The effect is small: the exit velocity is reduced 
only 1% with the off-optimum profiles. Because of this 
insensitivity to the p(x) profile, any reasonable approxi- 
mation to p(x) can be employed in the computer pro- 
grams. A convenient source of a p(x) profile for a given 
nozzle is the isentropic A(p) curve substituted into the 
A(x) profile for the nozzle. 

Figures 13, 14, and 15 illustrate some of the details of the 
flow in the 50-in. nozzle. The variation of gas and liquid 
velocity with distance is shown in Fig. 13. The initial 
velocities are 7.3 ft/s, but a slip velocity of 12 ft/s is 
established within the first inch of travel. The velocities 
then increase slowly to about 35 ft/s for the gas and 
20 ft/s for the liquid at a distance of 10 in. from the 
entrance. Rapid acceleration begins at this point, and 
the gas and liquid reach velocities of 235 and 165 ft/s, 
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Fig. 13. Variation of liquid and gas velocities 
with distance 

respectively, at the throat and 425 and 300 ft/s, respec- 
tively, at the exit. 

Figure 14 presents the temperature variation through 
the nozzle. Beginning with equal temperatures, the gas 
and liquid leave the nozzle with a 10°F temperature 
difference. 

The void fraction r a /{ 1 + r„) is shown in Fig. 15. The 
void fraction begins at 0.5, varies between 0.48 and 0.56 
during the initial stages of acceleration and drop breakup, 
and rises to 0.94 at the exit. 

V. Comparison With Experiment 

A. Nitrogen-Water Tests in the 50-in. Nozzle 

1, Construction and test procedure . An experimental 
nozzle was built whose contour was derived in Section IV 



Fig. 14. Variation of liquid and gas temperatures 
with distance 



Fig. 15. Variation of void fraction with distance 

(Fig. 6). Figure 16 is a drawing of the nozzle and injector, 
which were made of 6061 aluminum alloy, and Fig. 17 
shows the details of the injection tubes. The nozzle had 
an entrance diameter of 14.5 in. and converged at a 
20-deg half-angle to a diameter of 4.2 in. at a distance of 
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DIMENSIONS IN INCHES 


Fig. 17. Injection orifice details 


13.5 in. from the inlet. The nozzle then followed the 
tabulated curve to a throat diameter of 3.18 in. at 21.5 in. 
from the inlet, continued the tabulated curve to a diam- 
eter of 3.98 in. at 37.5 in. from the inlet, and diverged at 
a half-angle of 2.5 deg to an exit diameter of 5.2 in. at 
51.4 in. from the inlet. The nozzle was fabricated in three 
sections, which were bolted together and sealed with 
O-rings. A row of pressure taps was provided for measure- 
ment of the pressure profile. 

The injector consisted of 550 0.25-in. -diam tubes pro- 
jecting through slotted holes in the injector face. Water 
entered the injector through six 1.5-in. -diam inlets, flowed 
into the space between the middle and back plates, 
entered the injection tubes (Fig. 17), and left parallel to 
the injector face through three orifices at the tip of each 
tube. Nitrogen entered the injector through six 0.75-in.- 
diam inlets, flowed into the space surrounding the tubes 
between the middle plate and injector face, entered slots 
in the injector face, and flowed axially into the nozzle at 
three points around each tube in line with the water jets. 
The water orifices at the end of each tube consisted of 
slots 0.08 in. wide at three places around the circum- 
ference of each tube, with the upstream edge flush with 
the injector face; the slots in the initial test series were 
0.040 in. wide, with the upstream edge 0.040 in. from the 
injector face. The nitrogen orifices were three slots in 
the tube support holes, each 0.010 in. wide radially, 


occupying 60 deg of the circumference and extending 
1.0 in. axially to join with a 0.28-in.-diam entrance 
annulus. The slots in the initial test series were 0.005 in. 
wide. A 1.0- in. -diam bolt in the center extended through 
the three plates to reduce deflection. Pressure taps were 
provided in the water and nitrogen manifolds for deter- 
mining injection pressure. 

The nozzle was mounted on a parallelogram thrust 
stand, and the thrust was measured by a strain-gage force 
transducer mounted between the nozzle and a rigid 
support. Water and nitrogen flowed through tubes con- 
nected to overhead supply pipes through flexible hoses. 
Water was supplied from a pump capable of delivering 
180 lb/s of water at 500 psi. Nitrogen was supplied from 
twelve 25-ft\ 2200-psi cylinders at flow rates up to 10 lb/s. 
The flow leaving the nozzle was deflected downward into 
a sump from which the water returned to the pump while 
the nitrogen exhausted to the atmosphere. Figure 18 
shows the nozzle in operation during a test. 

The nozzle was operated at inlet pressures from 110 to 
230 psia, water flow rates from 90 to 180 lb/s, and nitro- 
gen flow rates from 2.5 to 9.5 Ib/s. The mean exit velocity 
V ranged from 250 to 540 ft/s, and the exit Mach number 
V(p(r«/p) ,/2 /( 1 + r a ) (Ref. 25, p. 7) ranged from 2.2 to 2.5. 
The nozzle sounded like an ordinary gas nozzle and had 
a 90-dB sound level. The pressure at the nozzle inlet 
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Fig. 18. Experimental nozzle in operation 


fluctuated by no more than 1 p si according to measure- 
ments made with a transducer having a frequency re- 
sponse to 1000 Hz. The exhaust jet was white and opaque, 
with some axial striations visible around the circum- 
ference. The jet had a sharp boundary as it left the nozzle, 
but the boundary became diffuse a few inches down- 
stream as a result of entrainment of atmospheric air. 

Figure 19 is a l-/*s flash photograph of the jet. The 
motion is stopped (0.003-in. travel at 300 ft/s), but the 
appearance differs little from that seen by the eye. 
The striations in the jet show more clearly; they are 
probably remnants of the injector pattern. 

The difference between the water injection pressure 
and the nozzle inlet pressure, the latter measured at a 
tap 3 in. downstream from the injector face, is plotted as a 
function of water flow rate in Fig. 20. The pressure drop 
was 14 psi at the maximum flow rate of 180 lb/s with the 


0.08-in. orifice width and 36 psi with the 0.04-in. orifice 
width, and the pressure drop was proportional to the 
1.5 power of the flow rate. The nitrogen injection pressure 
drop is shown in Fig. 21. The range of pressure drops for 
the tests was 10 to 60 psi for the 0.010-in. slots and 40 to | 

170 psi for the 0.005-in. slots. Liquid Freon 1301 (CBrF :t ) 
was used in place of nitrogen in one test, and Fig. 22 
shows the measured Freon injection pressure drops, 
which ranged from 7 to 25 psi. 

The flow rates given as rh\ and m y in Figs. 20-22 
should more strictly be labeled the component flow rates 
m b and m ay respectively. The solubility of nitrogen in 
water, however, is low enough that the differences are 
negligible, and it has been customary to employ rhi for m b 
and rh g for m u . Similarly, the mixture ratio r = rhi/rhy has 
customarily been employed where the component ratio 
r r = rhi/rha is really meant, and this practice will be 
followed in the remainder of the report. 
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Fig. 19. Exit jet at 1.0-fi.s exposure 
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NITROGEN INJECTION PRESSURE DROP A p f , psi WATER INJECTION PRESSURE DROP A/fc.psi 




Fig. 21. Nitrogen injection pressure drop 



The gas injection pressure drops and gas flow rates in 
the 50-in. nozzle were too low to provide atomization of 
a significant fraction of the water to drops of 0.02-in. 
diameter or less, thus making the nozzle dependent on 
drop breakup during expansion. Lack of atomization 
on injection was verified in an experiment in which a 
single injector element with 0.04-in. water orifices and 
0.005-in. nitrogen slots was operated over the range of 
nozzle mixture ratios inside a 3-in.-diam transparent tube 
maintained at 150-psi nitrogen pressure. The water jets 
were observed to be unaffected by the impinging nitrogen 
jets except for being deflected 5 to 10 deg away from the 
injector face. The flow from the injector of the 50-in. 
nozzle consisted, therefore, of a uniformly distributed, 
but only coarsely atomized, dispersion of water drops in 
nitrogen. A drop-flow, rather than bubble-flow, regime 
existed because the volume ratio of gas to liquid at the 
nozzle entrance ranged from 1 to 4. 

Based on the apparent lack of atomization at the inlet 
of this nozzle, the initial drop diameter appropriate for 
analyzing the nozzle with the real-flow program would 
be any value large enough to exceed We = 6 early in the 
nozzle and require drop breakup. Computer runs were 
made with various initial drop diameters, and a diameter 
of 0.05 in. was found sufficient for breakup in all cases, 
and for exit velocities that were independent of further 
increases in initial drop diameter, as illustrated in Fig. 8. 
The same approach was employed in selecting initial 
drop diameters for analyzing the other experimental 
nozzles to be discussed; a diameter was selected in each 
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case which gave results on the flat portion of the velocity- 
versus-diameter curve. 

2. Pressure profile. Pressures were measured at seven 
stations along the nozzle during runs made at 150-psia 
inlet pressure. The measurements made at a mixture ratio 
of 39.1 are compared with theory for r = 40 in Fig. 23, and 
show close agreement. The theoretical curve requires 
some additional explanation. The program does not give 
p(x) for a prescribed nozzle contour, but it does give 
area as a function of pressure, A(p), and the area, like 
the velocity, is insensitive to the p(x) profile supplied 
(see Fig. 12). Thus, the A(p) values printed are accurate 
for the actual nozzle, which has a p(x) profile differing 
little from the one prescribed in the program. A valid 
theoretical p(x) curve for the experimental nozzle can, 
therefore, be obtained by converting the printed A(p) to 
x(p) values, using the A(x) values for the experimental 
nozzle, and this procedure was the origin of the theoreti- 
cal curve in Fig. 23. The measured pressures are within 
1 psi of the theoretical curve. 



Fig. 23. Comparison of theoretical and 
experimental pressure profiles 


The theoretical and experimental variations of pressure 
with mixture ratio at a fixed station x = 29.46 in. are 
compared in Fig. 24. The measured values are within 
1 psi of theoretical from r = 15 to r = 70, falling below at 
low mixture ratios and above at high mixture ratios, 
probably in response to the nozzle exit behavior discussed 
next. 

Figure 25 shows the theoretical and experimental vari- 
ation of nozzle exit pressure with mixture ratio for inlet 
pressures of 150 and 230 psia. The experimental pressure 
was measured with a 0.04-in. -diam tap located 0.2 in. 
upstream of the nozzle exit. The theoretical exit pressure 



Fig. 24. Comparison of theoretical and experimental 
variation of pressure with mixture ratio 
21 in. upstream of exit 



Fig. 25. Comparison of theoretical and experimental 
variation of pressure with mixture ratio at exit 
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for 150-psia inlet pressure ranges from 16 psia at r — 10 
to 11 psia at r = 80, but the measured pressure is constant 
at 14.1 psia, the laboratory atmospheric pressure. At 230- 
psia inlet pressure, the theoretical exit pressure ranges 
from 22 psia at r = 10 to 14 psia at r = 80, but the 
measured pressure is constant at 16 psia up to the maxi- 
mum mixture ratio of 25 attainable at 230 psia. Thus, in 
practice, the exit pressure tends to remain equal to at- 
mospheric pressure, at least at the wall, even though the 
flow is nominally supersonic and independent of atmos- 
pheric pressure. 

3. Flow rate . Table 2 presents the theoretical total 
flow rates for isen tropic flow, real flow without friction, 
and real flow with friction for mixture ratios from 10 to 80 
at nozzle inlet pressures of 150, 190, and 230 psia. The 
flow rates are plotted as a function of mixture ratio for 
each pressure in Figs. 26, 27, and 28. 

The flow rate for real flow without friction is the 
highest; with friction, the flow rate is 2% lower, and for 



Fig. 26. Comparison of theoretical and experimental 
flow rates at 150-psia nozzle 
inlet pressure 


Table 2. Theoretical performance of 50-in. nozzle with water and nitrogen 1 


Inlet 


Isentropic- 

Real- 

Flow with 

Isentropic 

Real 

Velocity with 

pressure 

ratio 

flow rate 

flow rate 

friction 

velocity 

velocity 

friction 

Po 


m j 

mt 

m 6 

Vi 

V 


psia 


lbm/s 

lbm/s 

lbm/s 

ft/s 

ft/s 

ft/s 

150 

10 

78.5 

87.1 


642.6 

556.7 

539.8 


20 

106.0 

119.0 


476.2 

407.9 

393.4 


30 

126.0 

142.3 


400.3 

342.6 

329.5 


40 

142.0 

160.4 


355.1 

304.8 

292.4 


50 

155.5 

175.2 


324.5 

279.7 

267.9 


60 

167.1 

187.9 


302.2 

261.8 

250.2 


70 

177.2 

198.8 

194.6 

285.1 

248.2 

236.9 


80 

186.3 

208.6 

203.8 

271.5 

237.5 

226.4 

190 

10 

98.9 

108.4 

106.9 

676.1 

595.9 

575.8 


20 

132.8 

146.8 

144.3 

502.9 

439.8 

422.9 


30 

157.2 

174.2 

171.1 

424.2 

371.2 

355.8 


40 

176.5 

195.3 

191.7 

377.4 

331.4 

316.8 


50 

192.5 

212.7 

208.3 

345.8 

305.1 

291.0 


60 

206.2 

227.2 

222.5 

322.8 

286.2 

272.4 


70 

218.1 

239.9 

234.4 

305.3 

272.0 

258.4 


80 

228.7 

251.0 

245.2 

291.3 

260.8 

247.4 

230 

10 

1 19.1 

129.4 

127.5 

702.7 

627.3 

604.8 


20 

15?.l 

174.0 

171.0 

524.5 

465.7 

446.6 


30 

187.5 

205.2 


443.7 

394.7 

377.0 


40 

209.8 

228.9 

wmBmm 

395.8 

353.5 

336.7 


50 

228.1 

248.3 

242.8 

363.6 

326.3 

309.9 


60 

243.7 

265.0 

259.1 

340.2 

306.8 

290.8 


70 

257.2 

278.8 

272.5 

322.3 

292.1 

276.3 


80 

269.0 

291.2 

284.0 

308.2 

280.5 

264.9 

•T 0 = 60°f. p, > 

= 14.1 p*ia, A t — 7.942 in. 2 , = 0.05 

in. 
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Fig. 27. Comparison of theoretical and experimental 
flow rates at 190-psia nozzle 
inlet pressure 

isentropic flow it is about 10% lower than for frictionless 
real flow. The isentropic flow rate is lowest because, with 
no slip, the gas velocity is lowest and, with fixed throat 
area, the gas flow rate is lowest. At fixed mixture ratio, 
this results in the lowest total flow rate. 

Table 3 presents the experimental flow rates, and they 
are plotted with the theoretical curves in Figs. 26, 27, 
and 28. The experimental flow rates are about 3% lower 
than the theoretical values (real flow with friction) at the 
lowest mixture ratios tested and about 5% lower than 
theoretical at the highest ratios tested. The fact that the 
experimental flow rates are between the isentropic and 
real-flow values probably reflects the underestimation of 
liquid drag, discussed in Section II, because of the 
assumption of We = 6 diameter for all the drops and 
the use of the solid sphere drag coefficient. 

4 . Exit velocity. Table 2 also presents the theoretical 
exit velocities for isentropic flow, real flow without fric- 
tion, and real flow with friction. The velocities are plotted 
in Figs. 29, 30, and 31. The isentropic velocity is highest; 
the real velocity without friction is 10 to 15% lower, and 
friction reduces the velocity another 3%. 

The experimental exit velocities are given in Table 3 
and plotted with the theoretical curves in Figs. 29, 30, 
and 31. The experimental exit velocities are 3 to 5% 
higher than theoretical (real flow with friction), the differ- 
ence being attributed to the underestimation of drop drag. 



Fig. 28. Comparison of theoretical and experimental 
flow rates at 230-psia nozzle 
inlet pressure 

The experimental exit velocity is greater than theo- 
retical in this nozzle by an amount about equal to the 
velocity loss from friction. Thus, the real velocity without 
friction gives, fortuitously, the best agreement with 
experiment, predicting the measured values within 2%. 
Advantage was taken of this fact in the cycle studies 
presented in Ref. 27, where the nozzle exit velocities 
employed were the real-flow frictionless values. 

B. Freon— Water Experiment 

Six runs were made with liquid Freon 1301 (CBrF a ) in 
place of nitrogen to verify that the gas phase could be 
produced by contact vaporization at the inlet of a two- 
phase nozzle, an important requirement in MHD power 
systems. The results showed that complete vaporization 
occurred, in that the data agreed as well with the theory 
as in the nitrogen-water tests. 

1. Flow rate. Table 4 presents the theoretical and 
measured flow rates at the six mixture ratios and inlet 
temperatures obtained. The flow rates are plotted as a 
function of mixture ratio in Fig. 32. Real flow without 
friction is the highest, real flow with friction 2% lower, 
and isentropic flow another 4% lower. The experimental 
flow rates are about 2% lower than theoretical (real flow 
with friction). 

2. Exit velocity. Table 4 also presents the theoretical 
and measured exit velocities, which are plotted as a 
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Table 3. Measured performance of 50-in. nozzle with water and nitrogen 11 



Test 

series 

Inlet 

pressure 

P° 

psia 

Nitrogen 
flow rate 

m g 

Ibm/s 

Water 
flow rate 
m f 

Ibm/s 

Total 
flow rate 

Mixture 

ratio 

r 

Thrust 

F 

Ibf 

Exit 

velocity 

V 

ft/s 


1 

150 ± 5 


■HR 

98 

14.2 

■ 

468 






113 

19.4 

1 

405 





1 

128 

26.3 

1 

360 

I 

5 



4.04 


141 

33.9 


326 




3.49 

153 

156 

43.8 

1441 

297 

] 



2.73 

176 

179 

64.5 

1441 

259 




2.56 

184 

187 

71.9 

1464 

252 

I 

2 

150 ± 2 

6.03 

92 

98 

15.3 

1424 

468 

j 



5.75 

99 

105 

17.2 

1402 

430 




5.16 

109 

114 

21.1 

1448 

409 




5.07 

113 

118 

22.3 

1418 

387 




4.35 

123 

127 

28.3 

1418 

359 




4.32 

130 

134 

30.1 

1424 

342 




3.70 

145 

149 

39.2 

1424 

307 

- 



3.1 2 

161 

164 

51.6 

1424 

279 




2.71 

176 

179 

64.9 

1424 

256 

| 

3 

150 ± 5 

6.40 

90 

96 

14.1 

1450 

486 




4.65 

120 

125 

25.8 


371 




4.00 

135 

139 

33.8 

I 

336 




3.44 

150 

153 

43.6 


305 




2.97 

165 

168 

55.6 

1456 

279 

- 



2.83 

170 

173 

60.1 

1457 

271 


3 

190 ± 5 

7.45 

120 

127 

16.1 

1940 

491 




6.45 

135 

141 

20.9 

1916 

437 




5.60 

150 

156 

26.8 

1914 

395 




4.88 

165 

170 

33.8 

1915 

362 

- 

3 

230 ± 5 

9.40 

135 

144 


2425 

542 




8.40 

150 

158 


2425 

494 




7.39 

165 

172 


2420 

453 




7.05 

170 

177 


2440 

444 

— 










* T 0 = 60 ± 10°F, p 9 = 14.1 psia, A ( =7.942 In.* 


Table 4. Theoretical and measured performance of 50-in. nozzle 
with water and Freon 1301 at 150 psia® 


Inlet 


Isentropic- 

Real- 

Flow rate 

Isentropic 

KB 

Velocity 

Measured 

Measured 

temperature 

ratio 

flow rate 

flow rate 

with friction 

velocity 

I I 

with friction 

flow rate 

velocity 

To 

r 

m, 

m t 

m 5 

Vi 

n 


m t 

V 

°F 


Ibm/s 

Ibm/s 

Ibm/s 

ft/s 


ft/* 

Ibm/s 

ft/s 

46 

3.97 

■g 

MBM 

MEM 

406.4 

■a 

366.3 

mm 

372 

48 

4.66 

mbm 

mSSm 

mEM 

384.9 

MEM 

346.3 

■BP' 

355 

51 

6.38 

146.2 

mOm 

MEM 

344.8 

mSSM. 

309.5 

MEM 

315 

52 

7.1 1 

151.9 

160.7 

mEM 

332.0 

308.9 

297.8 

MEM 

314 

53 

8.79 


173.1 

170.0 

308.0 

286.9 

276.0 

169 

284 

54 

1 1.00 


187.0 

183.5 

284.9 

266.0 

255.3 

— 

— 

55 

13.66 

190.6 

26 l.l 

196.9 

264.6 

247.9 

237.4 

193 

247 

*p 0 = 150.0 ptfa, p e = 14.1 psia, = 7.942 in.* 

D 0 = 0.05 in. 
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Fig. 29. Comparison of theoretical and experimental 
exit velocities at 150-psia nozzle inlet pressure 



Fig. 31. Comparison of theoretical and experimental 
exit velocities at 230-psia nozzle inlet pressure 




Fig. 30. Comparison of theoretical and experimental 
exit velocities at 190-psia nozzle inlet pressure 


Fig. 32. Comparison of theoretical and experimental 
flow rates using Freon 


JPL TECHNICAL REPORT 32-987 


31 












values, showing that the loss due to friction is fortuitously 
counterbalanced by the drop drag underestimation. 

C. Six-Inch Nozzle 

Water-air experiments with a 6-in. -long two-phase 
nozzle at 514 psia inlet pressure and mixture ratios from 
4 to 33 were reported in Ref. 1. Figure 34 is a drawing of 
the nozzle, which had an inlet diameter of 1.62 in., throat 
diameter of 0.20 in., and exit diameter of 0.53 in. The 
two-phase, two-component program was employed to 
calculate the theoretical flow rate and exit velocity of this 
nozzle, using the isentropic pressure profile and 0.01-in. 
initial drop diameter, a diameter which caused breakup 
in the nozzle. 

L Flow rate . Figure 35 compares the theoretical and 
experimental flow rates. The experimental flow rates 
are 1 to 2% lower than with the theoretical values (real 
flow with friction). The fact that the flow rate is closer to 
theoretical than in the larger nozzle is attributed to the 
shorter time available for drop breakup, which was only 
a fraction of the drop natural period. The shorter exposure 
time would give a drop diameter distribution with a 
higher mean value, more closely approaching the Weber 
number = 6 value. 

2. Exit velocity . The theoretical and experimental exit 
velocities for the 6-in. nozzle are compared in Fig. 36 with 
the mean of the data represented by the dashed curve. 



Fig. 34. 6-in. nozzle and injector 


WATER AND FREON 130! 
150 psia 

THEORETICAL 

p o * 150 psia 
T 0 = 46-55 *F 
Oq * 0.05 in. 
p 0 s 14. 1 psia 
O MEASURED 

Pq - 150 ± psia 
T 0 * 46-55 °F 
p € * 14.1 psia 


8 300 [- 


X 

w 200 


^ — WITH FRICTION, V b 
■ ISENTROPIC, Vj | 


- REAL, V 


MIXTURE RATIO r 


Fig. 33. Comparison of theoretical and experimental 
exit velocities using Freon 


function of mixture ratio in Fig. 33. The exit velocity for 
real flow without friction is 7% below isentropic, and the 
velocity is reduced another 3% with friction. The mea- 
sured values are 3% above theoretical (real flow with 
friction) and again agree best with the real frictionless 
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EXIT VELOCITY, ft/s FLOW RATE . Ws 



MIXTURE RATIO r 


ig. 35. Comparison of theoretical and experimental 
flow rates for 6-in. nozzle 



MIXTURE RATIO r 


Fig. 36. Comparison of theoretical and experimental 
exit velocities for 6-in. nozzle 


The measured exit velocities range from 3% lower to 3% 
higher than with the theoretical values (real flow with 
friction). 

D. Steam Nozzle 

In Ref. 28, flow rates were reported for water-steam 
mixtures expanding through a converging-diverging 
nozzle, with the initial quality ranging from zero to 20% . 
The initial conditions and the measured pressure profile 
for the runs at 500-psia inlet pressure were supplied to the 
one-component program using an initial drop diameter 
of 0.01 in., which was large enough to cause breakup. 
Figure 37 compares the computed flow rates with the 
data from Ref. 28 (Fig. 4). The curve is the flow rate with 
friction; the flow rate without friction is 3% higher. 
The predicted flow rate is 30% below the measured value 
at a quality of 1.0%, and 15% below at 3% quality. 
Agreement is between 1 and 10% at higher qualities. 
The difference between predicted and measured flow 
rates at low quality is interpreted as the result of poor 
uniformity of the flow entering the experimental nozzle 
at those qualities. 



Fig. 37. Comparison of theoretical and experimental 
flow rates for steam-water nozzle of Ref. 28 
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E. Liquid Nitrogen Nozzle 

Experiments were reported by Richard (Ref. 29) in 
which room-temperature nitrogen gas was mixed with 
liquid nitrogen at the inlet of a two-phase nozzle. The 
injector provided a uniform spatial distribution of the 
liquid. The purpose of the experiments was to determine 
whether any gain in exit velocity over the value for equal 
gas and liquid temperature could be obtained with the 
initial temperature of the gas higher than that of the 
liquid. The measured velocities were such as to indicate 
that thermal equilibrium was reached on mixing, with 
no gain from the elevated temperature of the gas. 

To obtain nozzle inlet conditions for the one-component 
program, it was assumed that the injected gaseous and 
liquid nitrogen streams reached equal velocity and tem- 
perature immediately at the exit of the injector at the 
pressures measured there in the experiments; the pres- 
sures ranged from 206 to 289 psia. Under the assumption 
of thermal equilibrium, the resulting mixture consisted of 
subcooled liquid nitrogen in all cases. The velocity of the 
liquid, calculated from the flow rate and nozzle area at 
the injector exit, ranged from 104 to 114 ft/s. The liquid 
was then assumed to accelerate in the converging section 
of the nozzle until the saturation pressure, ranging from 
69 to 188 psia, was reached. The corresponding flow 
areas were all larger than the throat area, but the veloci- 
ties were such that the first increment of pressure drop 
below saturation would produce a percentage increase in 
velocity which was smaller than the percentage increase 
in volume from vaporization. Thus, the liquid could not 
vaporize upstream of the throat, and the only flow 
condition possible with thermal equilibrium was sub- 
cooled liquid upstream of the throat, saturated liquid at 
the throat, and vaporization starting immediately down- 
stream of the throat. The velocity at the throat could then 
be calculated from the flow rate, liquid density, and 
throat area. Using this velocity together with saturated 
liquid nitrogen properties as inlet conditions to the one- 
component nozzle program, the exit velocity was cal- 
culated, and the comparison with Richard's data is shown 
in Fig. 38^ The calculated exit velocity with friction 
included, Vs, agrees with the measurements within an 
average deviation of 4.3%. The good agreement is con- 
sidered to verify both the one-component program as 
applied to this case and the assumption of thermal 
equilibrium upstream of the throat in Richard's experi- 
ments. Of particular interest is the fact that friction must 
be included for good agreement and that the friction pro- 
duces a 14% velocity reduction, the largest of any of the 
examples considered. The cause is the larger wall area 



Fig. 38. Comparison of theoretical and experimental 
exit velocities for liquid nitrogen nozzle of Ref. 29 


per unit flow rate and lower average void fraction in 
the liquid nitrogen nozzle than in the others. Thus, the 
friction can account for the velocity reduction below 
isentropic, and it is not necessary to postulate a lag in 
vaporization as was done by Richard. The free-stream 
velocity V is only 3.5% below isentropic, because of the 
high density of the cold gas and the low surface tension 
of the liquid, resulting in good atomization at a low slip 
velocity. 

VI. Summary and Conclusions 

A one-dimensional analysis of two-phase nozzle flow of 
liquids and gases was made. The analysis employed five 
relations: the momentum and energy equations for the 
mixture, the drag and heat-transfer relations between the 
phases, and the Weber number — 6 criterion for drop 
diameter. The latter relation, and the use of the solid- 
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sphere drag coefficient, were the main simplifications and 
probably the main source of error in applying the analysis 
to experimental nozzles; in an actual nozzle, the drops 
have a range of diameters extending below the Weber 
number = 6 limit, and the drag coefficient can be higher 
than for a sphere because of flattening. 

Four computer programs were written for computation 
of real and isentropic flows of single-component and two- 
component mixtures. The real-flow programs provide the 
options of a prescribed pressure profile or an optimum 
pressure profile, and of selecting constant drop diameter 
or breakup to the Weber number = 6 limit. 

Two-phase nozzle design was investigated using the 
programs. It was found that the optimum nozzle profile 
has an elongated throat region, that the required shape 
is insensitive to drop diameter and nozzle length, and 
that the exit velocity is insensitive to departures from the 
optimum shape. For given fluids and pressures, the exit 
velocity was found to be insensitive to nozzle length and 
flow rate beyond certain minimum values; for water and 
nitrogen expanding from 150 to 14.1 psia at a mixture 
ratio of 40, for example, the exit velocity is between 
275 and 300 ft/s (78 to 85% of isentropic) for flow rates 
from 10 to 1000 lbm/s at the optimum lengths; for a flow 
rate of 160 lb/s, the exit velocity is between 275 and 
292 ft/s for all lengths between 5 and 225 in. This insen- 
sitivity of nozzle exit velocity to nozzle shape and size 
results from the inverse dependence of drop size on slip 
velocity, which apparently acts to adjust the slip to a 
value that holds the exit velocity in the range 75 to 90% 


of isentropic. Thus, it is difficult to design either a very 
good or a very poor two-phase nozzle. 

Theoretical and experimental values of flow rates and 
exit velocities were compared for two experimental 
nozzles, 50 and 6 in. long, operating with water and 
nitrogen, and water and air, respectively. For the large 
nozzle, the measured flow rates were 3 to 5% lower, and 
the measured exit velocities 3 to 5% higher, than the 
predictions. The deviations are attributed to the over- 
estimation of drop diameter and underestimation of drop 
drag coefficient in the analysis. For the small nozzle, the 
measured flow rates were 1 to 2% lower, and the mea- 
sured exit velocities from 3% lower to 3% higher, than 
the predictions. The closer agreement between measure- 
ment and prediction in the small nozzle is attributed to 
closer equality between the mean and the Weber num- 
ber = 6 drop size due to shorter time for drop breakup. 

Comparisons of theoretical and experimental results 
for one-component flow were made for two experiments 
reported in the literature. The theoretical flow rate for a 
steam nozzle differed from the experimental values by 
as much as 30%, but the experimental nozzle had no 
method of ensuring spatially uniform liquid distribution 
at the inlet, which is a requisite in the theory. The theo- 
retical exit velocity for a two-phase nitrogen nozzle 
agreed within 5% with the measurements. 

It is concluded that the analysis and computation 
method presented in this report will permit prediction of 
flow rates and exit velocities to an accuracy of 5% for 
two-phase nozzles having spatially uniform flow. 


Nomenclature 



Computer 


Computer 


Symbol 

Input or 
Output 
Name 

Definition 

0 7 7 Input or 

Symbol ' , , 

J Output 

Name 

Definition 

a 

— 

defined by Eq, (11) 

c ag CAG(P,T) 

specific heat of component A gas 

A 

A 

A 

AB 

flow area, in. 2 

flow area after velocity and tem- 


at pressure p and temperature T, 
Btu/lbm °R 



perature equalization, in. 2 

C<(i CAL(T) 

specific heat of component A liq- 

b 

— 

defined by Eq. (12) 


uid, Btu/lbm °R 
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Nomenclature (contd) 


Symbol 

Computer 
Input or 
Output 
Name 

Definition 

S 

CBG(P,T) 

specific heat of component B gas 
at pressure p and temperature T , 
Btu/Ibm °R 

Cb t 

CBL(T) 

specific heat of component B liq- 
uid, Btu/lbm °R 

Co 

— 

drop drag coefficient 

Ca m 

CDM 

value of C D at midpoint of a pres- 
sure step 

Cj 

— 

skin-friction coefficient for nozzle 
wall 

C/, 

CFIM 

skin-friction coefficient for inside 
wall of an annular nozzle 

c Sl 

— 

skin-friction coefficient for pure 
liquid flow 

C/ 0 

CFOM 

skin-friction coefficient for wall of 
a circular nozzle or for outside 
wall of an annular nozzle 

C 9 

— 

specific heat of gas, Btu/lbm °R 

Cg m 

CGM 

specific heat of gas at the midpoint 
of a pressure step, Btu/lbm °R 

Cl 

— 

specific heat of liquid, Btu/lbm 
°R 

Cl 

CL(P) 

specific heat of saturated liquid at 
pressure p (single-component 
flow), Btu/lbm °R 

C lm 

CLM 

specific heat of liquid at the mid- 
point of a pressure step, Btu/lbm 
°R 

D 

D 

drop diameter, in. 

^mai 

— 

maximum drop diameter, corre- 
sponding to We = 6, in. 

D 0 

DO 

initial drop diameter, in. 

F d 

— 

drag force on liquid, Ibf 

h 

HM 

heat-transfer coefficient between 
gas and drops at midpoint of a 
pressure step, Btu/h ft 2 


Symbol 

Computer 
Input or 
Output 

Name 

Definition 

K 

HG(P) 

enthalpy of saturated vapor at 
pressure p (single component 
flow), Btu/lbm 

H 

H 

inverse Henry’s Law constant, 
in. 2 /lbf 

K 

KAG(T) 

thermal conductivity of compo- 
nent A gas at temperature T , 
Btu/h ft °R 

K 

KBG(T) 

thermal conductivity of compo- 
nent B gas at temperature T, 
Btu/h ft °R 

h 

KGM 

thermal conductivity of gas phase 
at midpoint of a pressure step, 
Btu/h ft °R 

L 

— 

nozzle length, in. 

L 

L(P) 

latent heat of vaporization at pres- 
sure p (single-component flow), 
Btu/lbm 

L a 

LA(T) 

latent heat of vaporization of 
component A at temperature T, 
Btu/lbm 

Lb 

LB(T) 

latent heat of vaporization of 
component B at temperature T, 
Btu/lbm 

L m 

LM 

latent heat of vaporization at 
midpoint of a pressure step (single- 
component flow, Btu/lbm) 

m a 

— 

flow rate of component A, Ibm/s 

m« 

<> 

— 

flow rate of component A gas, 
lbm/s 

m ai 

— 

flow rate of dissolved component 
A } lbm/s 

mb 

— 

flow rate of component B, lbm/s 

rh bg 

— 

flow rate of component B gas, 
lbm/s 


— 

flow rate of dissolved component 
B, lbm/s 
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Nomenclature (contd) 


Symbol 

Computer 
Input or 
Output 
i Name 

Definition 

fn g 

MG 

gas flow rate, Ibm/s 

TTlg 

MGB 

gas flow rate after velocity and 
temperature equalization, Ibm/s 

rhi 

— 

isentropic total flow rate, lbm/s 

rhi 

ML 

liquid flow rate, lbm/s 

mi 

MLB 

liquid flow rate after velocity and 
temperature equalization, lbm/s 

rh t 

EMT 

total flow rate, lbm/s 

rhd 

— 

flow rate corrected for boundary 
layer at throat, lbm/s 

M 

— 

momentum flux, Ibm ft/s 2 

N 

— 

number flow rate of drops, s _1 

V 

P 

pressure, lbf/in. 2 

Pa 

PA 

partial pressure of component A, 
lbf/in. 2 

Pa 

PAB 

partial pressure of component A 
after velocity and temperature 
equalization, lbf/in. 2 

Pb 

PB 

partial pressure of component B, 
lbf/in. 2 

Pb 

PBB 

partial pressure of component B 
after velocity and temperature 
equalization, lbf/in. 2 

PK 

PBO(T) 

vapor pressure of pure component 
B at temperature T, lbf/in. 2 

Q 

Q 

separator parameter defined by 
Eq. (123) 

r 

R 

ratio of liquid mass flow rate to 
gas mass flow rate 

r 

RB 

value of r after velocity and tem- 
perature equalization 

r<x 

RA 

ratio of gas flow area to liquid 
flow area, piVi/rp (J V g 


Symbol 

Computer 
Input or 
Output 
Name 

Definition 

r c 

RC 

ratio of component B mass flow 
rate to component A mass flow 
rate 

r v 

RV 

ratio of gas volume flow rate to 
liquid volume flow rate 

n 

RVB 

ratio of gas volume flow rate to 
liquid volume flow rate after ve- 
locity and temperature equaliza- 
tion 

R 

— 

universal gas constant, 10.732 
(lbf/in. 2 )ftVlb-mole °R 

Re 

REM 

Reynolds number for drop motion 
relative to gas at midpoint of a 
pressure step, p g D V. s \J (i g 

Re f 

REF 

Reynolds number for liquid film 
flow at midpoint of a separator 

Ri 

— 

radius from axis to inner wall of 
an annular nozzle, in. 

Ro 

— 

radius from axis to outer wall of 
an annular nozzle, in. 

^0 

REDOM 

Reynolds number based on ve- 
locity boundary layer thickness at 
wall of a circular nozzle or on out- 
side wall of an annqlar nozzle at 
midpoint of a pressure step 


REDIM 

Reynolds number based on ve- 
locity boundary layer thickness on 
inside wall of an annular nozzle 

R$ 

— 

radius from axis to flow centerline 
of an annular nozzle, in. 


RAXO 

value of R<p at nozzle inlet, in. 

s 

S 

ratio of slip velocity to mean ve- 
locity V,/V 

So 

— 

optimum value of slip s 

s . 

SA 

Sutherland constant for compo- 
nent A, °R 

s„ 

SB 

Sutherland constant for compo- 


nent B, °R 
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Nomenclature (contd) 


Computer 
„ 7 7 Input or 

S « mbo1 Output 

Name 

Definition 

T 

T(P) 

saturation temperature at pressure 
p (single component flow), °R 

f 

TB 

temperature after velocity and 
temperature equalization, °R 

T. 

TG 

gas temperature, °R 


TGO 

temperature of gas at nozzle inlet, 
°R 

r, 

TL 

liquid temperature, °R 

Tio 

TLO 

temperature of liquid at nozzle 
inlet, °R 

T 0 

TO 

temperature at nozzle inlet (isen- 
tropic flow), °R 

i> ai 

— 

volume flow rate of component A 
liquid, ft 3 /s 

®»i 

— 

volume flow rate of component B 
liquid, ft 3 /s 

v„ 

— 

volume flow rate of gas, ft Vs 

V 

V 

velocity (isentropic flow), ft/s 

V 

VB 

mean velocity (V g + rVi)/( 1 + r), 
ft/s 

V g 

VG 

gas velocity, ft/s 


VGO 

gas velocity at nozzle inlet, ft/s 

Vi 

VL 

liquid velocity, ft/s 

Vlo 

VLO 

liquid velocity at nozzle inlet, ft/s 

v s 

VS 

slip velocity V g — V h ft/s 

v 6 

VBD 

mean velocity including boundary 
layer, ft/s 

Vo 

VO 

velocity at nozzle inlet (isentropic 
flow), ft/s 


WAG(P,T) 

effective molecular weight of com- 
ponent A gas at pressure p and 
temperature T 


WAGB 

effective molecular weight of com- 


ponent A gas after velocity and 
temperature equalization 



Computer 



c 7 7 Input or 
Symbol * 

* Output 

Definition 



Name 



w„, 

WAL 

molecular weight of component A 



WBG(P,T) 

effective molecular weight of com- 
ponent B gas at pressure p and 
temperature T 

j 

i 

W., 

WBGB 

effective molecular weight of com- 
ponent B gas after velocity and 
temperature equalization 

i 

i 

i 

i 

w bl 

WBL 

molecular weight of component B 


We 

— 

Weber number, p^V-D/Scr 


w s 

WG 

effective molecular weight of gas 


W g 

WGB 

effective molecular weight of the 



gas after velocity and temperature 
equalization 


X 

X 

axial distance from nozzle inlet, in. 

x, 

XG 

quality m g /m t 

X* 

XGB 

quality after velocity and temper- 
ature equalization 

X.0 

XGO 

quality at inlet of nozzle 

\)i 

YI 

distance from flow centerline to 
inner wall of an annular nozzle, in. 

Vo 

YO 

radius of circular nozzle, or dis- 
tance from flow centerline to outer 
wall of an annular nozzle, in. 

a 

ALPHA 

mass fraction of component A dis- 
solved in liquid 

a 

ALPHB 

mass fraction of component A dis- 
solved in liquid after velocity and 
temperature equalization 

p. 

BETA 

mass fraction of component B va- 
por in gas 

p 

BETAB 

mass fraction of component B va- 
por in gas after velocity and 
temperature equalization 

8 

— 

velocity thickness of boundary 
layer, in. 


f= - - 


38 


JPL TECHNICAL REPORT 32-987 


Nomenclature (contd) 



Computer 


Computer 

„ , 7 Input or 

S « mb ° l Output 

Definition 

„ 7 7 Input or 

S « mb ° l Output 


Name 


Name 

8* 

— 

displacement thickness of the 
boundary layer, in. 

A r — 

A s — 

Si 

DELI 

velocity thickness of boundary 
layer on the inner wall of an an- 
nular nozzle, in. 

i 

t-r 

< 

s? 

DELSI 

displacement boundary layer 
thickness on inside wall of an an- 

A T g — 



nular nozzle, in. 


S m T 

— 

difference between T g and T* at 
midpoint of a pressure step, °R 

A T, — 

So 

DELO 

velocity thickness of boundary 
layer on wall of a circular nozzle 
or on outer wall of an annular 

AT, — 



nozzle, in. 




1 

N 

> 

<d 

K 

DELSO 

displacement boundary layer 
thickness on wall of a circular 
nozzle or on outer wall of an an- 
nular nozzle, in. 

AV; - 

S l T 

— 

difference between T g and Ti at 
beginning of a pressure step, °R 

A Vf - 

8,V 2 

— 

difference between V 2 and V'\ at 
beginning of a pressure step, ft 2 /s 2 

Ax — 

ST 

— 

difference between T g and T i 

9 — 

AA 


change in flow area across a pres- 
sure step, in. 2 

9i THI 

AH* 

— 

change in enthalpy across a pres- 
sure step for process n, Btu/lbm 
°R 

9 io THIO 

A "S 

— 

change in flow rate of component 
A gas across a pressure step, lbm/s 

0 o THO 

A S 

— 

change in flow rate of component 
B gas across a pressure step, lbm/s 


Am„ 

— 

change in gas flow rate across a 
pressure step, lbm/s 

K THOO 

Ap 

DP 

pressure change in one pressure 
step, lbf/in. 2 



Definition 


change in r across a pressure step 

change in slip across a pressure 
step 

change in gas temperature across 
a pressure step, °R 

change in gas temperature due to 
velocity and temperature equaliza- 
tion, °R 

change in liquid temperature 
across a pressure step, °R 

change in liquid temperature due 
to velocity and temperature equal- 
ization, °R 

change in square of mean velocity 
across a pressure step, ft 2 /s 2 

change in V 2 across a pressure step, 
ft 2 /s 2 

change in V‘f across a pressure step, 
ft 2 / s 2 

change in X across a pressure step, 
in. 

momentum thickness of boundary 
layer, in. 

momentum boundary layer thick- 
ness on inside wall of an annular 
nozzle, in. 

initial momentum boundary layer 
thickness on inside wall of an an- 
nular nozzle, in. 

momentum boundary layer thick- 
ness on wall of a circular nozzle or 
on outer wall of an annular nozzle, 
in. 

initial momentum boundary layer 
thickness on wall of a circular 
nozzle or on outer wall of an an- 
nular nozzle, in. 
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Nomenclature (contd) 


Computer 


Symbol 

Input or 
Output 

Name 

Definition 

A 

ALAM 

Lagrangian multiplier for nozzle 
contour optimization, ft/s 2 

/S 

VIAG(T) 

viscosity of component A gas at 
temperature T, Ibm/ft h 

P<*i 

VIAL(T) 

viscosity of component A liquid at 
temperature T, lbm/ft h 


VIBG(T) 

viscosity of component B gas at 
temperature T, lbm/ft h 

Pb t 

VTBL(T) 

viscosity of component B liquid at 
temperature T, lbm/ft h 

Pff 

VIG(P) 

viscosity of saturated vapor at 
pressure p (single-component 
flow), lbm/ft h 

Pg 

VIGM 

viscosity of gas at midpoint of pres- 
sure step, lbm/ft h 

Pi 

VIL(P) 

viscosity of saturated liquid at 
pressure p (single-component 
flow), lbm/ft h 

PI 

VILM 

viscosity of liquid at midpoint of a 
pressure step, lbm/ft h 

Pi 

VILB 

viscosity of liquid after velocity 
and temperature equalization, 
lbm/ft h 

P a l 

ROAL(T) 

density of component A liquid at 
temperature T, lbm/ft 3 

' Pbj 

ROBL(T) 

density of component B liquid at 
temperature T, lbm/ft 3 

Pg 

ROG 

density of gas, lbm/ft 3 

Pg 

ROG(P) 

density of saturated vapor (single- 
component flow), lbm/ft 3 

Po 

ROGB 

density of gas after velocity and 
temperature equalization, lbm/ft 3 

Ps m 

ROGM 

density of gas at midpoint of a 
pressure step, lbm/ft 3 

Pi 

ROL 

density of liquid, lbm/ft 3 


Symbol 

Computer 
Input or 
Output 

Name 

Definition 

Pi 

ROL(P) 

density of saturated liquid (single- 
component flow), lbm/ft 3 

Pi 

ROLB 

density of liquid after velocity and 
temperature equalization, lbm/ft 3 

Phn 

ROLM 

density of liquid at midpoint of a 
pressure step, Ibm/ft 3 

a 

SIGMA 

surface tension of liquid, dyn/cm 

<7 

SIG(P) 

surface tension of saturated liquid 
at pressure p (single-component 
flow), dyn/cm 

( 7 

SIG(T) 

surface tension of component B 
liquid at temperature T, dyn/cm 

Tw i 

TWIM 

shear on inner wall of an annular 
nozzle, lbf/in. 2 


TWOM 

shear on wall of a circular nozzle 
or on outer wall of an annular 
nozzle, lbf/in. 2 


PHI 

angle between axis and flow cen- 
terline of an annular nozzle, deg 

IDi 

WIM 

angle of inner wall relative to flow 
centerline for an annular nozzle, 
deg 

ti) 0 

WOM 

angle of wall relative to axis for a 
circular nozzle, or angle of outer 
wall relative to flow centerline for 
an annular nozzle, deg 

Subscripts 


e 


nozzle exit 

i 


isen tropic flow 

m 


midpoint of a pressure step 

t 


nozzle throat 

0 


nozzle inlet 

1 


beginning of a pressure step 

2 


end of a pressure step 
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